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1 .  Abstract 


This  report  describes  work  done  during  the  period  5/1/86  - 
4/30/87  on  three  areas  of  MPD  research:  (a)  A  rigorous  theoretical 
examination  of  the  one-dimensional  flow  of  a  self-field  accelerated 
plasma,  clarifying  the  roles  of  area  change,  sonic  vs.  magnetoacoustic  choking 
and  finite  magnetic  Reynolds  number.  (b)  Continued  development  of 
computat ional  codes  for  axisymmetric  MPD  flows,  and  ( c )  Results  of  a 
first  test  series  on  two  channels  designed  to  verify  predictions  on 
the  effects  of  area  variation,  showing  ability  to  redistribute  current 
by  this  means. 
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( a3 )  Extension  of  the  Above  to  Applied  Field  MPD  Flows . 

Confined  to  axisymmetric  conditions,  with  no  current  spokes. 

(bl)  Experimental  Study  of  Various  Trhuster  Configurations.  The  aim  here 
was  to  validate  theoretical  concepts  developed  previously  about  the  effects  of 
electrode  contouring. 

(b2)  Diagnostics  Development.  To  apply  previously  acquired  optical  and 
signed  processing  equipment  for  the  measurement  of  various  plasma  parameters  in 
connection  with  (bl)  above. 

In  the  following  sections  we  give  an  overview  of  the  accomplishments  in 
these  tasks  and  of  the  occasional  departures  from  the  original  plan.  Reference 
will  be  made  to  three  papers  presented  by  our  group  at  the  1987  International 
EP  Conference,  where  the  details  of  most  of  the  work  were  described.  These 
papers  are  included  here  as  Appendices. 

2 .  Outline  of  Work  Accomplished 

2.1  Theory  In  our  efforts  to  clarify  the  dynamics  of  self- field  accelerated 
MPD  flows,  we  had  previously  developed  the  limiting  concept  of  a  high  Magnetic 
Reynolds  Number  (Rm)  plasma  (l),  which  was  useful  in  understanding  the  close 
connection  between  flow  area  variation  and  distribution  of  current  along  the 


length  of  an  MPD  thruster  of  slender  geometry,  and  also  made  clear  the  dominant 
role  of  magnetoacoustic  waves  and  magnetoacoustic  choking  phenomena. 

The  practical  range  of  Rm  is  from  2  to  10  or  jc- ,  which  is  not  extremely 
high,  and  it  was  not  entirely  clear  to  what  extent  these  "high  Rm"  results 
would  apply  to  realistic  conditions.  Particularly  troublesome  to  some 
colleagues  was  the  fact  that  in  this  limit,  thermal  effects  were  relegated  to  a 
definitely  secondary  role,  with,  for  instance,  no  mention  of  the  delicate 
balance  required  at  a  sonic  point,  as  had  been  discussed  by  several  previous 
authors  (2),  (3).  In  fact,  these  sonic  conditions  were  thought  to  be  entirely 
responsible  for  the  determination  of  the  all-important  channel  voltage  for  a 

given  current  (2),  and,  following  this,  for  the  occurrence  of  conditions,  at 
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high  values  of  (current)  /(mass  flow  rate),  where  the  back  emf  exceeded  locally 
the  channel  voltage,  a  condition  that  was  identified  as  leading  to  the  symptoms 
known  as  "onset"  (3). 

We  were  also  aware,  as  was  pointed  out  in  Ref.  (1),  that  the 
approximations  leading  to  the  simple  behavior  of  the  "high  Rm"  plasma  would 
break  down  in  regions  of  high  current  concentration,  near  the  inlet  and  exit  of 
a  channel.  However,  the  proper  way  to  account  for  these  regions  was  not  known 
prior  to  the  work  reported  here.  Even  at  this  time,  the  complications  brought 
about  by  the  simultaneous  effects  of  high  magnetic  Reynolds  number  and  high 
Hall  parameter  in  two  or  three-dimensional  geometries  are  not  well  understood. 
For  example,  each  of  these  two  effects  contributes  to  the  existence  of  a 
current  concentration  at  the  root  of  the  cathode,  while  the  Hall  effect 
alleviates  the  corresponding  high  Rm  concentration  at  the  anode  root.  While 
these  factors  are  separately  understood  now,  no  simple  theory  of  the  combined 
effect  exists  yet.  In  view  of  the  difficulties  mentioned  above,  which  were 
severe  even  with  no  Hall  currents  (as  in  a  one-dimensional  channel  with 
segmented  electrodes  in  a  Faraday  connection),  we  decided  to  postpone  work  on 
multi-dimensional  formulations  and  concentrate  our  efforts  on  a  precise 
analysis  of  the  one-dimensional  case.  The  results  of  this  work  reported  in 
Ref.  (4).  which  is  included  here  as  Appendix  1. 
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The  key  result  is  that  even  as  Rm-)”  ,  the  effects  of  magnetic  diffusion 
(i.e.,  ohmic  dissipation)  are  to  be  retained  in  ever  thinner  inlet  and  exit 
layers,  in  order  to  satisfy  the  boundary  conditions  on  the  magnetic  field  (for 
the  prescribed  current).  This  is  entirely  analogous  to  the  classical  very  thin 
viscous  boundary  layers  in  nearly  inviscid  fluids.  And,  also  in  analogy  to 
boundary  layer  theory,  the  effects  of  such  layers  may  be  important  for  some 
purposes  (as  for  drag  in  the  viscous  case)  while  they  may  be  ignored  for  others 
(as  for  lift).  In  the  MPD  case,  we  show  in  Appendix  1  how  the  existence  of 


these  sharp  inlet  and  exit  current  layers,  where  dissipation  concentrates,  is 
of  minor  consequence  for  voltage  prediction,  but  strongly  affects  calculations 
of  dissipation  (hence  efficiency  and  thrust).  It  is  also  shown  that  sonic 


passsage  does  occur  within  the  inlet  layer,  and  that  proper  satisfaction  of  the 
smooth  sonic  passage  conditions  is  indeed  essential  if  this  layer  is  to  be 
resolved  properly.  However,  the  ma.jor  flow  metering  function  is  still  found  to 
be  fulfilled  by  magneto-sonic  choking  at  a  channel  throat,  well  downstream  of 
the  sonic  point. 

In  addition  to  clarifying  these  points  and  hence  for  the  first  time 


providinng  a  corrct  solution  for  high  Rm,  the  work  reported  in  App.  1  explored 
numerically  the  range  of  validity  of  the  Rm->”solut.ion.  It  was  found  that,  for 
channels  that  are  contoured  so  as  to  evenly  distribute  the  current  (which  leads 
to  increased  efficiency  and  reduced  electrode  stress),  the  high  Rm  results  hold 
well  for  Rm  down  to  about  3  or  so,  and  depart  only  gradually  at  lower  Rm.  In 
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contrast,  in  constant-area  channels,  the  strong  dissipation  near  the  exit  at 


finite  Rm  produces  positive  pressure  gradients  that  slow  the  flow  and  reduce 
the  local  magnetic  convection,  leading  to  thermal  choking  at  Rm-5,  to  embedded 
shocks  between  Rm=5  and  Rm-2,  and  to  subsonic  flow  at  lower  Rm.  These  regimes 
had  been  previously  reported  by  King  (2).  Thus,  high  Rm  predictions  for 
constant-area  cases  are  of  limited  validity  for  Rm  <8.  On  the  other  hand, 
constant  area  channels  are  probably  impractical  due  to  the  strong  current 
concentrations  and  the  attendant  reduction  in  thrust  efficiency. 

Finally,  we  compared  the  voltage  and  thrust  predictions  of  our  models 
(with  finite  Rm)  to  those  reported  by  the  Princeton  group  (5)  for  an 
experimental  slender  channel,  and  found  good  agreement  in  both  cases,  including 
a  good  account  of  the  transition  between  the  electrothermal  and  electromagnetic 
regimes . 

2.2  Numerical  Simulation.  Preliminary  results  of  a  developmental  axisymmetric 
MPD  flow  numerical  model  were  included  in  our  previous  yearly  report  (6).  Thus 
work  has  been  continued  and  refined,  with  the  most  recent  results  being 
discussed  in  Ref.  7  and  in  our  Appendix  2. 

The  refinements  have  been  motivated  by  numerical  difficulties  encountered 
at  the  higher  current  levels,  for  a  given  flow  rate.  They  have  included  a 
smoother,  numerically  generated  curvilinear  grid  (see  Fig.  1  of  App.  2)  with 
concentrations  near  corners  and  walls,  as  well  as  a  more  rigorously  calculated 
Jacobian  matrix,  including  both,  flow  and  electromagnetic  equations.  The 
Jacobian  matrix  is  the  matrix  of  partial  derivatives  of  the  various  governing 
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equations  w.r.t.  the  chosen  dependent  variables,  evaluated  at  all  the  grid 
points,  and  is  used,  after  inversion,  to  generate  a  new  linearized  solution  to 
all  the  equationns.  These  are  re-linearized  and  the  process  is  iterated.  It 
has  been  found  that  this  Newton-Raphson  method  converges  very  rapidly  (4-6 
iterations,  compared  to  thousands  of  time  steps  in  most  time-marching  schemes1, 
but  requires  large  amounts  of  computer  storage.  This  means  that  the  grid 
resolution  must  necessarily  be  limited,  for  any  given  machine  capacity.  We  are 
currently  working  with  20  x  50  grids  in  our  Micro-Vax  computer.  This 
relatively  coarse  grid,  coupled  with  the  very  steep  property  gradients  found  in 
certain  regions  of  the  solution,  leads  to  the  occurrence  of  numerical 
fluctuations  which  necessitate  an  objectionable  amount  of  artificial  numerical 
damping,  and  which,  at  high  currents,  induce  negative  density  "pockets", 
followed  by  generalized  breakdown  of  the  solution.  It  must  be  also  stated  that 
the  stability  limits  of  the  Newton-Raphson  algorithm  for  the  kind  of  mixed 
convective-diffusive  and  elliptic-hyperbolic  type  of  equations  we  are  using  are 
not  known  at  present,  and  hence  other  factors  may  be  at  work  in  limiting  the 
method. 

Within  these  limitations,  the  results  obtained  for  6  g/sec  of  Argon  at 
currents  to  10  KA  are  quite  interesting  (see  Appendix  2).  The  strong  plasma 
density  concentration  at  the  cathode  tip  (plasma  focus)  is  clearly  indicated 
in  Fig.  6  of  App.  2,  as  is  the  existence  of  rarefied  regions  near  the  anode  lip 
and  at  the  cathode  shoulder.  The  cathode  tip  and  root  current  concentrations, 
well  known  to  experimenters,  are  visible  in  Fig.  3  of  App.  2. 

Incidentally,  the  numerical  instability  appears  to  be  unrelated,  at  least 

2 

in  principle,  to  the  expected  anode  depletion  at  high  I  “/m,  since  the  plasma 
density  over  most  of  the  cathode  is  only  some  20%  lower  than  that  at  the 
cathode  at  I  =  10  KA,  and  is  well  above  that  level  required  to  supply  the  local 
random  electron  flux  to  the  anode  wall. 


-7- 


$ 

ft 

& 

»  **• 
l\*. 


si 

$ 

* 


& 


R 


Work  is  underway  now  to  exploit  the  larger  capabilities  of  a 
multiprocessor  vectorized  computer,  of  which  partial  ownership  was  acquired  by 
a  special  AFOSR  supplementary  grant.  Two  lines  of  attack  are  planned:  (a)  to 
explore  the  behavior  of  the  Newton-Raphson  code  with  finer  grids  and  (b)  to 
utilize  sophisticated  time-marching  techniques  (Flux-Corrected  Transport 
Codes),  now  available  in  our  VAX.  The  latter  work  will  be  undertaken  by  a  new 
doctoral  student,  Mr.  Eli  Niewood  who  is  supported  by  an  NSF  fellowship,  while 
the  Newton-Raphson  work  will  be  continued  by  Mr.  J.  Marc  Chanty. 

The  slower  than  expected  development  of  our  numerical  work  has  forced  us 
to  postpone  consideration  of  applied-field  cases.  This  is  still  contemplated 
for  the  future,  since  no  essential  new  difficulties  are  expected  once  the  self¬ 
field  case  is  better  in  hand  (with  the  exception  of  spoke  generation,  which 
would  require  some  form  of  3-D  calculation) . 

2.3  Experimental  Work.  In  cooperation  with  Dr.  Peter  Turchi  and  collaborators 
of  the  Washington  office  of  R  &  D  Associates,  we  have  designed,  built  and 
tested  several  MPD  channels  in  order  to  learn  about  the  degree  of  fidelity  of 
the  necessarily  simplified  theories  we  have  developed  previously.  This  work 
was  largely  carried  out  by  doctoral  candidate  D.  Heimerdinger ,  who  spent 
several  months  at  the  R&D  lab.  in  the  course  of  the  development  of  the  required 
experimental  apparatus  and  facilities.  He  was  also  helped  by  a  second  graduate 
student,  Mr.  D.  Kilfoyle  wo  installed  the  optical  diagnositcs  and  participated 
in  the  experiments  themselves. 

The  preliminary  results  of  tests  with  two  channels,  one  of  constant  area 
and  one  featuring  a  divergent  section,  are  reported  in  Ref.  8  which  is  also 
reproduced  in  Appendix  3.  Additional  results  with  a  convergent-divergent 
channel  will  be  reported  later,  together  with  a  more  complete  analysis  of  the 
data  presented  here. 


a/Iv ./ -\.V  .xvV-V wxvmVh".  jATV’jV..  yvx.V,. 


seas^MaiiieaiBiiteassasasaaStssaes^  agawas  asssaa  j&SSSE&SSm, 


-8- 


The  tests  are  pulsed  (  0.5  msec),  allowing  a  relatively  simple 
apparatus.  In  an  attempt  to  avoid  the  complicating  efforts  of  axial  symmetry, 
the  channels  were  built  with  a  fairly  large  cathode  radius  and  small  cathode- 
anode  gap,  thus  approaching  two  dimensionality.  This  also  has  the  effect  of 
boosting  the  level  of  current  that  can  be  passed  before  "onset".  On  the  other 
hand,  it  also  leads  to  weakened  electric  and  magnetic  fields  for  a  given  total 
current,  and,  on  the  initial  test  series,  some  difficulty  was  encountered  to 
obtain  smooth,  symmetric  discharges  at  the  design  conditions. 

What  apeared  to  be  symmetric  discharge  conditions  (as  determined  by  sub¬ 
microsecond  frontally  taken  photographs)  were  finally  obtained  by  a  combination 
of  increasing  the  current,  decreasing  the  flow  rate,  and  simply  accumulating 
over  one  hundred  pulses  in  order  to  condition  the  electrode  surfaces  by  burning 
away  the  inevitable  adsorbed  and  oxydized  layers.  The  conditions  obtained  this 
way,  which  are  those  reported  in  App.  3,  were  well  above  the  theoretically 
expected  onset,  but  during  the  test  series  with  the  constant  area  channel  we 
did  not  encounter  the  MHz  voltage  fluctuations  which  have  come  to  be  regarded 
as  the  signature  of  onset.  Upon  examination  of  the  channel  at  the  end  of  the 
tests,  it  was  clear  that  large  arcs  had  been  present  on  both  anode  and  cathode, 
and  they  were  visible  in  some  of  the  high-speed  photographs  taken  during  the 
tests.  It  was  only  later  that  we  have  been  able  to  show  conclusively  that 
onset  was  indeed  exceeded  in  these  tests,  and  this  has  come  mainly  from 
monitoring  the  near-anode  voltage  drops  in  a  subsequent  test  series  (results  to 
be  reported  later).  The  absence  of  voltage  "hash"  in  the  constant -area  channel 
which  initially  misled  us  into  believing  we  were  still  below  onset,  has  been 
seen  to  be  a  peculiarity  of  this  geometry,  despite  other  clear  evidence  of 
anode  depletion.  The  contoured  channels  do  exhibit  voltage  "hash",  although 
with  some  delay  and  with  individually  differentiated  patterns. 
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Given  the  fact  that  conditions  during  the  tests  reported  here  were  above 


onset,  we  cannot  expect  to  find  good  correlation  with  theoretical  results  that 
are  for  smooth,  non-depleted  operation.  There  are,  however,  clear  indications 
of  the  main  effect  that  was  pursued,  namely,  a  significant  redistribution  of 
current  due  to  the  exit  flaring  in  the  second  channel,  in  the  sense  of  reducing 
the  strong  exit  concentration  measured  in  the  constant-area  channel  (see  Fig. 

13  of  App.  3) . 

As  indicated  before,  diagnostics  development  has  proceeded  in  parallel 
with  the  rest  of  the  experimental  program,  and  spatially  and  spectrally 
resolved  data  of  the  near-exit  plume  were  obtained  during  some  of  the  tests. 
Additional  measurements  were  made  using  an  R&D  plasma  generator.  These  results 
are  not  given  here  in  detail,  however,  since  their  interpretation  is 
incomplete  at  the  reporting  time.  We  will  simply  note  that  we  use  a  fast 
digitizing  Silicon  Intensified  Target  (SIT)  camera,  together  with  a  1.26  m 
Czerny-Turner  spectroscope,  to  obtain  at  once  complete  transverse  profiles  of 
the  plane,  with  the  spectrum  either  resolved  over  a  region  containing  several  H 
lines  for  Te  measurements,  or  containing  the  resolved  HS  line  for  Stark 
measurement  of  n^  .  Measurement  of  T  j  by  Doppler  broadening  is  also 
attempted.  An  algorithm  for  automatic  Abel  inversion  of  the  data  has  been 
developed,  which  uses  FFT  spectral  methods,  allowing  pre- f i l ter i ng  of  the 
generally  noisy  data. 

3.  Personnel  Associated  with  this  Research 

The  work  was  directed  by  Associate  Professor  Manuel  Mart inez-Sanchez ,  with 
the  cooperation  of  Graduate  Research  Assistants  D.  Heimerdinger  (Doctoral 
Candidate),  J.  March  Chanty  (Doctoral  Candidate.,  D.  Kilfoyle  M.S.  Candidate' 
and  Eric  Sheppard  (Doctoral  Candidate  . 
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THE  STRUCTURE  OF  SELF-FIELD  ACCELERATED  PLASMA  FLOWS 
By  Manuel  Martlnez-Sanche2 
M.I.T.,  Dept.  of  Aeronautics  h  Astronautics 


Abstract 

An  analysis  is  presented  of  self-field 
acclereted  que»s  i-one-d  imens  1  one  I  piasirnna  f  low* 
wi»h  zero  axial  current.  '.‘he  dominance  of 
magnetic  convection  over  diffusion  leads  *o  the 
formation  of  inlet  and  exit  conncenf rat  ion  layers 
which  are  analyzed  by  asymptotic  matching.  The 
sonic  passage  condition  is  found  to  control 
injector  conditions,  but  not  channel  voltage. 
Throats  are  shown  to  be  approximately  magneto- 
snnic.  Performance  factors  accounting  ‘or  inlet 
convergence  and  nozzle  divergence  are  derived  bv 
exploiting  t be  smallness  of  pressure  forces.  The 
effect  of  finite  magnetic  Hevnolds  nuumbers  is 
seen  to  be  moderate  down  to  'bn  -  3  for  ronntoured 
channels,  but  to  Rm  -  10  for  cons fan t-ar°a 
chnnne'.s.  Contouring  is  shown  to  be  effective  for 
control  of  current  and  dissipation  concentrations. 

1.  !nt rorluction  This  work  is  nn  attempt  to 
clarify  the  roles  of  various  fluid  and 
electromagnetic  effects  on  determining  the 
performance  and  the  distributions  of  current  end 
density  in  M^O  self-field  thrusters.  The  general 
problem,  with  ail  its  atomic  physics  implication 
and  with  due  account  for  real  electrode  effects, 
is  clearly  beyond  simple  analysis,  and.  must  be 
tackled  by  a  combination  of  experiment  and 
numerical  simulation.  However,  even  the  idealized 
flow  problem,  where  only  minimal  physics  is 
included,  has  so  far  eluded  analysis  to  a  large 
extent.  Since  even  at  the  simplest  meaningful 
level  of  abstraction  (i.e.,  quasi-one  dimensional 
channel  flow)  the  fluid  mechanics  of  a  r“!f-fieid 
accelerated  plasma  is  complex  and.  unfamiliar,  it 
is  worthwhile  to  go  through  its  analysis  tn  ns 
comolete  a  fashion  aa  possible.  The  results  will 
shed  light  on  several  aspects  of  MPU  operation, 
including  the  role  of  thermal  choking  to 
determining  thruster  voltage',  the  possibility  of 
the  back  e.m.f.  exceeding  ;he  applied  voltage"  the 
role  of  channel  contouring  and  the  usefulness  of 
the  large  Magnetic  Reynolds  number  approximation,^ 
as  well  as  the  low  gas  pressure  approximation, 
both  of  which  could  be  exploited  in  multi¬ 
dimensional  modeling  work. 

?.  formulation  We  wish  to  postulate  conditions 
under  which  a  quasi-one  dimensional  flow  model  is 
•enable.  for  a  low-pressure  partially  magnetized, 
flow  ns  in  self-field  wpq,  this  is  complicated  by 
•he  Hall  effect,  which,  'or  ecjuipotent  iai  wails, 
produces  substantial  axial  currents  and  subsaquent 
’•-anxverse  density  gradients,  leading  eventually 
to  deplation  of  the  anodic  region.  ,,J* 

*This  work  was  supported  by  the  Air  force  Office 
nf  Scientific  Research  under  Grant  No.  AfOSH-HS- 
OllHO 

Associate  Professor,  Dept.  Aeronautics  and 
As  * ronaul ics , 

Member  A1AA 


Conceptually,  at  least  two  approaches  can  be  taken 
to  bypass  this  problem:  (a)  assume  a  high  enough 
power  level  that  the  Hall  effect  becomes 
negligible,  or  '  b !  assume  an  appropriate  axial 
electrical  field  can  be  imposed  through 
segmentation,  such  as  to  cancel  the  axial  'Hall) 
current.  Option  la)  exists  because  the  Hall 
parameter,  d,  can  be  shown  to  scale  ns  t  ,  trr 
for  a  given  fluid  and  geometry,  while  Baksht's 
•heory  indicates  that,  at  the  initial  (onset)  con¬ 
dition  is  u  l  1l"  ,  or,  assuming  operation  near 
onset,  c'U’/y 7  .  Option  'b'  allows  uniform 

•rnnsverse  conditions  even  in  the  presence  of  a 
finite  Hall  effect,  and  will  be  assumed  in  this 
work.  The  required  f  axial  field  can  be  found  a- 
posteriori,  and  no  attempt  will  be  made  here  to 
discuss  implementation  of  the  implied  segmented 
construction,  other  than  to  point  out  that  the 
fields  involved  are  typically  of  the  order  of  b-20 
V/  cm  and  should  pose  no  axial  breakdown  problems. 

The  steady-state  quasl-one  dimensional 
equations  of  mass  and  inviscid  momentum 
conservation,  written  m  terms  of  cross-section 
aversged  properties,  ere 


- 0  (2) 

where  9  is  the  magnitude  of  the  self-induced 
magnetic  * ie Id ,  which  is  perpendiculsr  to  the 
centerline  direction.  i'he  total  enthalpy 
conservation  law  is 


pu^-(h-*v>  -  ?.*T> 
f  cx  2 


where  the  average  current  density  J  is,  by 
asaumot ion,  purely  transverse  and  is  related  to 
the  magnetic  field  by 


The  condition  Vx?=0  can  be  shown  to  reduce 
in  our  case  of  zero  axial  current  to 

TE  H)  -  0  ;  F.  H  •  V  -  const.  (5) 
dx  v  v  y 

which  \*  the  sum*  aa  would  apply  in  the  absence  of 
Hall  effect.  Here,  M  •_  *  ♦hr'  \  n*.  er -a  lac*:  rode 
channel  heigh*:.  There  a,  however,  a  non-zero 
axiel  field  given  {on  average'  by 

(6) 

p  «  B  (F.  -uR) 

x  v 

where  8  m  the  Hell  parameter. 


yrujer  theae  condition*,  the  mean  electrical 
work  v..  j  become® 
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and  tha  energy  aquation  (Sq.  (3)),  combined  with 
Bqe.  .5)  and  (7),  integrates  to 

const. 


A(h  ♦  t-)  —  Vy3 

u« 


(8) 


where  W  la  tha  channal  width,  which  in  aaauaad 
constant.  Note  that  this  is  also  the  name  fora 
that  would  apply  with  zero  Hail  effect. 


To  complete  the  foraulatton,  the  transverse 
component  of  Ohm's  law  (neglecting  electron 
diffusion  effects'  is 


0  (E 


ufl ) 


(9) 


where  c  is  the  scalar  conduct iv ity .  Also,  an 
equation  of  state  is  needed,  most  usefully  in  the 
fora  h  =  h.P.o).  In  this  work  we  will  use  either 
a  simple  perfect  gas  model,  with  constant  ratio 
of  specific  heats,  or  an  ionization  equilibrium 
model  for  a  noble  gas. 


In  our  analysts  we  will  take  the  point  of 
view  of  having  a  proscribed  total  current  t 
(implying  a  given  80 =  Ugl/W)  and  flow  rate  d, 
with  a  geometrical  throat  of  area  A*.  A!,so  .given 


is  the  inlet  total  enthalpy, 


The  voltage  Vy  (or  the  field  E  at  some  section) 
is  then  to  be  determined.  Also  unknown  are  the 
inlet  density  and  velocity  (but  not  their  product, 


since  ai  =o,u0Ao  is  known).  As  to  exit 


conditions,  there  are  two  logical  alternatives: 

(a)  Since  the  exit  velocity  ib  likely  to  be  higher 
than  the  apeed  of  propagation  of  flow  disturbances 
(i.e.,  as  we  will  see,  the  magnetoacoust lcs 
speed',  we  can  leave  exit  conditions  unspecified, 
m  the  understanding  that  complete  expansion  to 
zero  magnetic  field  and  pressure  will  occur 
external  to  the  channel.  (b)  In  the  context  of  a 
qua»i-one  dimensional  model,  we  can 
bottle  up"  this  external  expansion  by  imposing, 
somewhat  artificially,  zero  magnetic  field  at  the 
exit  plane.  This  procedure  will  increase  the 
concentration  of  current  in  the  exit  region,  and 
will  thus  lead  to  stronger  Ohmic  losses  and  an 
underestimation  of  efficiency,  but  it  has  the 
advantage  of  providing  simple  and  consistent 
est ’metes  of  overall  performance  without  tha  need 
for  two  or  three-dimensional  external  flow 
calculations. 


It  is  convenient  to  no  non-dimensionalize 
variables  using  reference  values  derived  from  the 
prescribed  quantities.  These  reference  quantities 
are: 
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An  important  non-dimensional  parameter  is  the 
Magnetic  Reynolds  number,  Rn  -  juju  .,1. 

t  ret 

It*  inverse. 


Is  numerically  in. the  range  0. 1-0.3  for  realistic 
Magneto  Rlnsmadynamic  (MPU)  thrusters. 


In  what  follows,  w«  will  work  with  thm  non- 
dimensional  variables 


b  -  S/B,  , 

u  ■ 

u/u  ,  p 

raf 

’  p/pr.f 

P  -  D/or„ 

.  H 

•  h/h  ,  , 

ref 

l  •  t  H  , 
y  ref 

J  -  j/«  , 

*  "  ref 

,  a 

-  A/A*  -  H/H* 

,  K  •  x/L 

V  -  V  /V  „ 

y  ref 

(12) 

In  te rms  o f 

these 

variables,  the 

governing 

equations  become 

cua  ■  1 

(13) 

--  Hu 

c  u  — 
if. 

* 

b2)  -  0 

(  14) 
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v1 .  d>; 

2 

2€J  ; 

h-h(p,  Q) 

(15) 

cJ  -  S 

-  ub 

(  16) 

J  -  - 

dh 

c\r 

(17) 

Ea  - 

V 

(  18) 

Notice  that  Rq.  (15)  is  the  internal  energy 
equation,  obtainable  by  subtracting  the  product  of 
u  times  Rq.  (2)  fro*  the  differential  form  of  Kq. 
(8).  Clearly,  the  alternative  is  the  non- 
dimensional  form  of  Rq.  (8)  itself: 


-2 

+  V-  +  2Vb 


hto 

-2 


2  V 


(19) 


la 

In 


where  -  r  , 

V  "  h  -r  u 
to  o  2  o 

Since  J  -  -db/df|,  the  small  parameter 
seen  to  multiply  magnetic  field  derivatives, 
the  limit  as  e->  o  the  system  ( 13)  —  ( IB)  is 
therefore  singular,  with  possible  layers  of 
rapidly  varying  properties  at  the  channel’s  inlet 
and  exit,  due  to  the  presence  of  prescribed 
boundnrv  conditions  there,  and  with  the  body  of 
the  c.nanne',  between  sucn  layers,  constituting 
what  mathematically  speaking  would  be  called  the 
"outer"  region,  in  which  the  c.  terms  are 
insignificant. 


3 .  Bin gulnr  perturha tion  Problem  -  The 


Zero' th  Order _ Solution 


rn  the  limit  when  r 


the  width 
and  !  •  1 
ength  of 
to 


of  the  smgu'ar  "inner"  layers  near  f  =o 
fends  to  zer-o,  and  *  or  nearly  the  whole 
"he  chiinne ,  Kqa .  (15'  ano  (  .6)  reduce 


'ft 

_3 


(20) 

(21) 


Here  the  superscript  ( 1  )  refers  to  the  "outer" 
solution.  Ho.  ■  HCf 1  indicates  lsentroplo  flow  in 
"us  '-eg i0n .  Kq.  '?!),  combined  with  5a  m  y 
and  ?’u  a  =  1 , 

-  (22) 


gives 

cons  t . 
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"here  V  is  the  '  non-dlmenslonal)  channel  voltage 
fo  order  rero.  Notice  ►hat  this  voltage  is  the 
seme  throughout  t  ha  channel,  Including  the 
singular  "inner”  layers,  hut  may  Include  higher 
order  perturbation  tanas  for  non-sero  £. 


i 

iCK 
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The  nature  of  the  zero’ bb  order  "outer" 
solution  can  be  readily  seen  by  defining  an  • 
equivalent  pressure  which  indiudea  both  <na  and 
magnetic  contributions 


In  terms  of  P  .  the  equation  of  motion  (to  any 
order)  is  todd-ed"  a  0,  end  since  \f\  this  zero' th 
order  out«r  solution  both  p  find  b  have  been 
found  depend  on  ,>  «lnn«,  the  flow  la  barotropic 
n nd  behaves  the  same  way  ns  an  ordinary  Ideal  ga* 
flow  with  a  special  pressure-dens i ty 
relationship.  In  MPO^ flows  b  >>  p,  so  that,  to  a 
first  npproxiraat ion, w  ,  and  the  effective  ”t" 

of  the  flow  is  about  2.  The  speed  of  propagation 
of  smalt  disturbances  in  such  a  flow  is  (in  non- 
d ifens tone  l  terms)  _ _ _ 

3  /TO,’  (24) 

where  c  is  the  ordinary  speed  of  sound,  and 
2V  =  2  ~  can  be  identified  ns  the  square  of 

(be  A’.fven  velocity,  and  is,  in  fact,  the  dominant 
contribution  to  Eq.  (2b).  It  is  a  known  property 
of  ideal  gas  flow  that  the  gns  velocity  at  a 
*hront  must  equal  the  sma 1 1 -disturbance 
prooagation  velocity,  provided  the  downstream 
pressure  is  low  enough.  Therefore,  we  deduce  the 
existence  of  a  magnetosonic  flow  choking 
condition,  which,  for  a  prescribed  flow  rate  and 
throat  area,  must  imply  a  certain  relationship 
among  the  upstream  conditions.  In  our  problem 

"upstream"  conditions  must  be  those  immediately 
pest  the  thin  inlet  leyer,  through  which  we  can 
■live  nearly  discontinuous  variations.  !,et  us 
denote  >hi«  location  by  a  aubscrlpt  f.o 

distinguish  it  Jrom  inlet  conditions  (subscript 
'o','.  Since  "O.can  be  found  f;rom  the  continuity 
equation  and  FT;  depends  on  ft  !  and  o.  ,  the 
remaining  independent  variables  at  (1)  are  (-c  ,  b  , 
and  uj  In  addition,  the  voltage  V  is  still  an 
unknown.  The  condition  that  the  throat  be 
magnetosonic  is  then  of  the  form  *'(p  j,b  :,u  f  ;  V,' 

-  0,  the  specific  details  depending  on  trie  fluid 
model  adopted. 

0  0  _0  — 

Additional  equations  for  p  ,b  ,u  and  V  can 
be  found  by  asymptotic  matching  of ’th<*  "outer" 
solution  so  tar  discussed  to  the  "inner"  solution 
'hrcuvh  the  inlet  leyer.  Since  the  thickness  of 
*his  aycr  is  of  order  j,  we  can,  in  the  zero’ th 
approximation,  ignore  area  variation  in  it;  on  the 
orher  hand,  the  t  terms  in  Eqs.  (15)  and  0.6)  are 
now  to  be  retained.  It  is  convenient  >;o  re-scale 
distance  by  defining  an  inner  coordinate 


-  c/t  (25) 

The  "Inner"  equations  of  momentum, 
total  energy  and  Ampere’s  law  follow  from  (14), 
:  i'J!  and  '  lb)  respectively  ; 

1  du1  ^  d  (p*  +  b1  )  -  0 


h1  +•  u1*  +  2vbl  -  ii  +  2V 

: -  CO 


dbl  -  -  (V  -  QS1)  ;  e1  -  V_ 


Initial  conditions  for  the  equations  can  be 
applied  et  inlet,  as  has  been  done  in  (27)  (b^to) 

„  1  ,  ( hl  ♦  uli/;)*h  The  inlet  velocity 

u Vo'  ■  ’T,  anc*  preseWe  p  Ho)  *  p(  are  not  a- 
priori  known,  however. 

Eq.  (26)  integrates  to 
ul“  u0  +  a  j  p‘)  +  1-b12’  (29) 

Two  equations  for  th*  condition*  at  ^  M  can 
be  *imply  * ound  bv  expressing  the  Integrated  inner 
law*  1  l'-  )  and  (2‘i)  at  their  outer  limit: 


u,  +  a  0 1  p  0  -  p  *  i  -  b^  ' 


l/2(G1V+  2Vb  *  -  ht0  +  27 


An  additional  matching  condition  can  be 
obtained  from  Eq.  (2H)-  since  b'  must  smoothly 
tend  to  b:  as  fhe  inner  variable  r d 
tends  to  inf  inity,  we  must  have 

V  -  a.S’bJ  (32' 

This  is  in  fact,  equivalent  to  Eq.  (21),  when 
expressed  at  the  initial  point  of  the  outer  range. 

It  can  be  seen  now  that  Eqs.  (30),  (31)  and 
(32),  plus  the  magnetosonic  choking  condition,, 
will  determine  properties  at  ( L)  and  voltage  V 
provided  the  miet  velocity  ia  known  (the  inlet 
pressure  p  would,  then  roiiow  from  the  given  c'aU0 
and  Thl*  indicate*  that  some  additional 

condition  remain*  to  be  found,  and  thl*  will  bo 
seen  to  be  the  smooth  passage  through  the  sonic 
point  (which  occur*  wall  upstream  o*‘  the 
magnetos pn^ic  point,  in  fact,  within  the  inlet 
Laver'.  The  importance  of  this  condition  has  been 
emphasized  by  ^ing~  and  lawless  and 
Subramanian  **'  who  state  that  it  is  the  sonic 
passage  condi t ton  that  governs  the  channel 
voltage.  On  physical  grounds,  this  would  seen 
unlikely  (rr  fact,  Nef.  (3)  speaks  of  a  "non- 
mtuitive  result”'.  Indeed,  the  voltage  is  largely 
the  induced  back  e.m.f.,  determined  bv  the 
eventual  ga*  velo«  ify,  wh i  •  h  Is  many  * ime*  larger 
than  f he  speed  of  sound,  and  one  would  expect 
downstream  cpndittons  to  play  the  major  role  in 
determining  V.  That  this  is  ►he  case  can  be  seen 
from  the  discussion  so  rsr  •  the  initial  velocity 
u$  is  s  very  smal 1  quantity  (of  order  0.0’,  ns  our 
results  will  shew'  nnd  its  direct  contribution  to 
Eqs.  '  .:0 '  and  '3l)_can  be  safely  ignored.  Chi  the 
other  hand,  since  -  l,  setting  u„  =  o 

would  imply  an  infinite  density,  and  hence 
pressure  p^;  »n  ’act,  the  sma ’ l ,  but  non-zero  u* 
leads  to  a  rather  large  ,  but  •  U '.  1  to  a  smal  l 
pm  1  of  order  0.  !  .'  a*  nr  responds  *o  the  physical 
regime  of  operation  of  M},3  devices  (gas  pressure 
<<  magnetic  pressure'.  It  follows  that  both  u, 
and  p3  could  be  ignored  'and  hence  also  '  m 
Eqs.  39'  and  31  ,  and  V  could  be  calculated  to  a 
rai r  approximation  with  no  reference  to  nny  gas 
thermodynamics,  and  no  need  ror  a  sonic  pasaage 
condition.  As  we  will  explicitly  show  in  a  simple 
case,  the  role  of  this  sonic  passage  condition  is 
mostly  to  determine  the  precise  va*ue  of  ►he  inlet 
variables  ,  Tig,, etc.)  once  the  voltage  V  t* 

known  fro*  the  downs tren*  conditions. 


Th®  for*  of  th#  sonic  passage  condition  cm 
ba  seen  for  instance  by  solving  the  set  of 
equations  (13)  to  (17)  for  the  individual 
derivatives,  this  will  in  fact  ba  useful  for 
later  numerical  computations.  For  the  velocity 
derivative,  we  obtain 


_2  _2 

(u  -c  )  du 
*  d\ 


2  (?.  -  Zb)  E-dJb 


where  c  =/  V?  is  the  normalized  isentroptc 
speed  of  sound,  And 


The  existence  of  an  inlet  discontinuity  in  the 
magnetic  field  distribution  for  a  constant-area 
channel  (and  hence  of  an  inlet  current 
concentration)  when  R->  -»  was  already  suggested  in 
Ref.  4,  but  Its  magnitude  was  erroneously  given  as 
l-b,  =  1-2/3  —  i/3.  The  error  arose  froa  the 

identification  of  u,  with  u  ,.  In  the  same  paper, 
the  inlet  discontinuity  was  ignored  when  dealing 
with  a  convergent-divergent  channel.  wa  will  see 
that  this  is  justified  onlv  In  the  limit  of  very 
large  inlet  contraction  ratio,  although  It  may  be 
a  reasonable  approximation  for  practical 
contraction  ratios. 


K 


(34) 


is  s  thermodynamic  quantity  which  reduces  to 
in  the  case  of  an  ideal  gas  with  constant  specific 
heat  rstio  f  It  is  clear  that  when  j  =  c  ,  the 
whole  rig^t  hand  side  of  F.q.  '  113 >  must  be  s 
simultaneously  zero.  In  the  present  zero' th  order 


npprox imat ion  (t<'ll,  only  the  second  term  is  of 
concern  (which  is  equivalent  to  i gnor i ng  area 
variation  in  the  " inner" ^region ) .  The  ’wo  roots 
of  thia  second  tern  are  1  =  ub  and  P.  =  Kub,  and 
the  first  of  these  is  reached  at  the  outer  edge  of 
Che  initial  layer.  Since  ub  Increases  fron  near 
zero,  and  K>1,  the  condition 


E  *  (Ku'bS 

S  S 


(35) 


is  reached  inside  the  layer,  and  must  be  the  sonic 
passage  condition.  Here  the  subscript  a  indicates 
the  sonic  point.  Together  with  G'  =  c  and  with 
the  conservation  equations  expressed  a?  the  sonic 
point  (Eq.  (27',  (29)),  taia  allows  elimination  of 
all  variables  at  s,  yielding  the  final  required 
equation  for  the  calculation  of  all  properties  at 
stations  I.O)  and  (1)  and  of  tho  voltage 


4.1  tamp! e_7ero^_t h_0rder_9esu It s  He  hnva 

1  mp  emeu  t  ed  ‘Vi  nl».v«  pr  O'  e, lures  i.tr  Toth,  a 

oone»Miil-A  gas  m,,. tel  and  an  equilibrium  Argon 
model  al'.iaving  first  ionization  only,  and  verified 
*heae  modela  against  each  othor  in  low  tempera' ure 
eases.  rwo  examples  will  be  presented  here,  one 
for  a  convergen* -divergent  channel  with  e  =A  ,/A  -2 
and  the  other  for  a  constant -area  channel  (a  =  1). 
!n  both  cases  we  assumed  m/A  •  0.5  Kg/sec, 

R  =0. 1  Tesla  ,  T  =  400K. 

toe  ,o 

Thia  means  reference  quantities  uref=  7960  m/sec, 

pref=  •l‘,H0  •  nref~-  6-2fl  * 

2  . -  79b  V/m,  and  also  h  =  0.0032H. 

ret  to 

Table  1  ahows  the  principal  results  for  the 
convergent-divergent  channel.  Several  features  of 
these  results  are  noteworthy: 

a)  The  flow  is  sub-megnetosonic  (but  mostly 
supersonic)  in  the  "chamber”  pert  of  the  channel. 

b)  There  is  a  week  current  concentration  at  inlet 
(carrying  1  -  0.9401,  or  (A,  of  tne  total 
current) . 


Once  these  values  ere  known,  channel  exit 
conditions  must  be  calculated  by  solving  e  second 
set  of  "inner'  equations  similar  to  Kqs .  (26)-(28) 
for  the  rapid  variations  that  occur  near  the 
channel  exit,  where  the  condition  b  -  a  i.s  to  be 
imposed.  The  task  is  now  simple,  since  the 
•lgenvslues  of  the  problem,  i.e.,  V  and  u0  are  by 
his  time  known  quantities;  physically,  if  the 
exit  speed  is  greater  than  (or  equal  to)  the 
magnezoaonlc  speed,  exit  conditions  cannot 
inf luenca  upstream  conditions.  He  thus  can 
calculate  exit  velocity,  and  hence  thrust, 
e'ficiency  and  specific  impulse  of  the  'hruster  at 
the  specified  current  and  mass  flee*. 

The  thrust,  normalized  by  P  As ,  is 

ref 

r  •  +  Pe*.  (35a) 

where  the  second  tern  contains  no  magnetic 
contribution  because  of  our  assumption  that  8=0  at 
and  beyond  the  exit  plane.  "he  torus  ter 
efficiency  is  defined  here  as 

"  '  Mlntmijm_pqwer_.  for  a  given  thrust  and  flow 

thr  Actual  power 


V' 


or,  in  terms  of  th#  dimensionless  variables. 


i’V 


c)  Because  of  this  last  feature,  the  actual  Argon 
"emperature  In  the  chamber  is  only  3100  K,  too  low 
for  equilibrium  ionization.  Von-equt 1 ibrium 
effects  should  be  accounted  for,  but  this  result 
emphasizes  the  "off-loading"  of  the  channel  inlet 
area. 

d)  There  is  a  drastic  density  drop  through  even 
this  weak  current  layer.  This  wtll  hamper  efforts 
st  straight  numerical  simulation. 

e)  The  ratio  p  of  pressure  to  msgnetic  pressure 
remains  at  «  few  percent  or  below  throughout  the 
flow. 

for  comparison,  Tab’.e  2  shows  the  results  for 
a  non-contoured,  cons* *mt -area  channel  with 
otherwise  the  same  data.  The  entire  channel 
except  *  or  the  inlet  and  exit  "layers")  is  now 
Fmgnetosonic.  The  inlet  region  carries  (  1  - 
9. ”3',  or  J7\  of  *he  current,  the  rest  being 
carried  by  the  exit  current  concentration.  The 
extra  inlet  dissipation  raises  the  Argon  channel 
iemprrnture  to  7f;b/i  g  ionization,  =  1.18) 

Vi'  it  and  tha  exit  concentration1  lowers  the 
efficiency,  as  sb.sxn  -n  'he  ex  i  '  nnm'  imens  iona  l 
'Trust  of  about  1.  .6  1  from  'M  effects,  0.  16  from 

e 1 ec* ro' henna  1  expansion',  compared  to  !.5\  for 
‘he  -as"  in  "ab'e  1.  'n  fsct,  for  the  perfect 
expansion  of  "sb'e  ',  *  he  theemal  energy  resulting 
f'om  ‘he  Inlet  dissipation  is  fully  recovered,  and 
’he  thrust  efficiency  is  unity,  hut  even  if  we 
assume  no  recovery  beyond  ’he  throat,  an 
e f f i '  ency  of  0.95  is  calculated  'ignoring,  of 
-nurse,  all  transport  losses'.  By  comparison,  the 


'  35b  * 


thrust  efficiency  of  the  constent-area  channel  is 
only  0.730.  Not*  also  the  even  wore  extreme 
density  drop  through  the  inlet  layer  in  this  case. 

Notice  that  in  both  example*  above  the 
channel  length  is  not  specified,  beyond  stating 
that  it  euat  be  larg®  enough  to  sake  the  Magnetic 
Reynolds  nuffber  large.  For  Q  =  4CCJ014  */■  we  have 
fusing  Eq.  (11)),  R^-  40L,  with  l  In  #.  A  .’0  cs. 
Tong  ►hrustcr  would  thus  have  Rm  *  8,  or  = 

0.k’5.  We  will  see  in  Sec. 5  that  this  value  is 
*mnll  enough  to  ensure  a  reasonably  good 
appro**  im«  t  ion  from  our  zero’  order  solution, 
including  all  the  main  qualitative  features 
discussed. 

4.2  E xplicit  Zcro'th  Order  Results 
Neglecting  Pressure  .  It  was  stated  tn  Sec.  3  f ha t 
Hfw  pressure  forces  play  a  secondary  role  in  MF9 
thrusters  < except  as  discussed  in  Section  5.1,  tor 
low  Wm,  constant  area  cases)  and  this  is  to  some 
extent  verified  by  the  low  values  of  p  :n  the 
results  presented.  ft  is  therefore  of  interest  to 
simplify  the  formulation  by  dropping  p  from  the 
momentum  equation  (Eq.  14).  This  allows  us  to 
’gnore  the  internal  energy  equation  (Eq.  \5), 
whose  role  is  to  determine  p  consistently,  and 
also  to  drop  any  consideration  of  sonic  passage, 
since  the  absence  of  nnv  thermal  effects 
essentially  implies  a  very  large  Inlet  Mach 
number.  We  could  now  leave  the  inlet  velocity  u 
as  an  arbitrary  parameter,  or,  provided  we  0 

renounce  calculation  of  p^,  simply  set  it  to 
zero.  To  illustrate  the  smallness  of  the  effect 
of  a  finite  uq  ,  we  will  carry  it  ns  a  parameter  in 
the  calculations. 

Neglecting  ?>,  Eq.  (14)  reads 
-  du  b  db  n 

U  ~~T  +2  -  "0  /  or\ 

dt.  7  di.  (36) 

and  in  the  "outer"  region,  where  =  V,  this 
integrates  to  c 


uo“  <(u,)  ■*,4V(b,-b°)  (37) 

At  'he  same  tine,  V  =  u°  D  e.  so  that  the  magnetic 
*ield  b*  in  this  outer  region  is  related  to  the 
channel  area  through 


l/bJ  -i-4?(b’-b5  |  , 

Imposing  that  a  =  1  when  da/db  =  0  (throat) 
is  equivalent  to  impoatng  magnetosonic  choking. 
This  gives  'he  condition  , 

;  9  i/l  - 1  .  -  V 

^  *^b,  j'-*/TV  •  Ul'~T  (39) 

In  the  "inner"  region,  equation  (29)  reduces 
now  to  _ (  it 

u  •  u  *■  a  (1-b  ) 

o  o 


/  (x+ll 

*»*  '  l-x) 

and  than  tha  voltage  la 

5  -  1  x*l  * 


and  the  velocity  at  station  (1)  ia  given  by 

U  1  a,  U  -  *)  (4<i> 

An  » 

Assigning  vslues  to  x  between  1  and  3  we  thus 
obtain  t hm  results  shown  in  Table  3,  covering  the 
rengm  from  no  inimt  contraction  (Aj  =  A*)  to 
mf’.ni'  e  contraction  1 A^  > 'A  ).  For  a^  =  1  and  s,, 

=  2  these  rmsults  can  be  compared  to  those  in 
Tables  1  and  2,  which  include  the  pressure 
forcms.  For  a,  7  1,  neglecting  p  gives  V  ..  0.393 
if  we_take  uo  =  0.0147  Tnble  2),  or  V  =  0.285 
with  uo  -  0;  the  result  in  Table  l  is  7  =  0.158. 
Het ter  approximation  is  obtained  for  a-j  =  2  *v  = 
9.54J  Using  u0  =  0.014J,  versus  0.586  from  Table 
n.  It  seems  clear  that  neglecting  thermo  1 
effects  does  give  the  correct  trends,  particularly 
when  we  compare  other  quantities  os  well. 

To  calculate  the  exit  conditions  we  need  to 
specify  the  exit  area  ratio  af  =  \  /A  .  The 
magnetic  field  Just  prior  to  the  exit  current 
concent  rat  ion  is  not  yet  zero;  denoting  this  0 
location  by  the  subscript  (2)  we  can  calculate  b2 
by  substityting  a  =  a  into  Kq.  (38)  and  then 
u,  =  v/a  b-  •  Integrating  the  momentum  equation, 
in  the  4  ora  of  Eq.  (36),  across  the  exit  "inner" 
layer  then  gives  for  the  exit  velocity  (b  =  0): 

~Je  '  “2  +  ae(bi)  (46) 

Tt  is  noteworthy  that  this  exit  velocity  can 
be  explicitly  written  down  as  the  product  of  two 
factors,  one  deoending  upon  contraction  ratio,  a0 
-  ,  the  i  -  her  deoending  upon  expansion  rntio 

ae  r  Ae/A*.  rhe  result,  after  some  algebra,  ts 

.  t  ra0'?^(ae)  (47) 


Y  4-  1 

C  m  (-■ - 

0  ' 

i-i  x 


(x  given  by  Rq.  43)  (48) 


*  .  ^ - -  A  (TT  4.  2  COS  *  )  (49) 

e  o  a 

e 

The  efficiency  of  the  thruster  reduces  In 
this  limit  to 

n  1/2  mG2  G  2 

— vT“*  ‘  -50’ 

Combining  Eqs.  '47)  and  (44),  this  can  also 
be  expressed  in  'he  present  approximation)  as  the 
product  of  ‘ha  factors 


which,  when  substituted  into  the  matching 
condition  (32)  gives 

a~  '  bi  5’A,+a,(l-bJ  )!  '41) 

Eliminatiryf  7  between  (31)  and  (41)  gives  an 
equation  ror  b  ns  a  function  of  a,.  The  solution 
ia  best  expressed  Inversely.  Defining  an 
intermediate  variable  x  through 


b.y: 


n  "  n  (a  'a  ^ 

o  o  e  e 


e  .) 


■(^  ): 


Tha  factor*  f  and  n'  are  thu*  the  exit  velocity 
and  efficiency  of  a  purely  convergent  channel  with 
a  variable  aiaount  of  pre-throat  convergence,  while 
and  neare  the  velocity  and  efficiency  of  a 
purely  divergent  channel  of  varying  expansion  area 
ratio.  The  limiting  values  (uaing  u  =  o)  are 
shown  in  Table  4.  ° 


Notice  the  strong  beneficial  effect  of  Inlet 
convergence:  divergence  poet  'he  throit  is  also 
beneficial,  but  to  a  lesser  extent.  Even  though 
no  energy  equation  hna  been  used,  the  implication 
of  an  efficiency  less  than  unity  is  that  some  of 
the  electrical  energy  IV,  has  been  lost  as  thermal 
energy  in  the  exhaust.  The  efficiency  limit  of  1 
in  the  strongly  convergen t -d i vergen t  esse  is  due 
to  the  oliaination,  in  that  '.  isi’,  of  'be  s’rong 
current  concentrations  both  af  inlet  and  at  exit; 
it  is  these  dissipative  layers  'hat  are 
responsible  for  the  mef f icienc tea  in  the  other 


A  word  of  caution  is  in  order  here.  We  have 
ignored  the  physical  processes  Leading  to 
ionization  at  the  channel  inlet.  The  nmount  of 
inlet  contraction  must  be  carefully  balanced 
between  the  desire  1  or  high  efficiency  and  *  he 
need  for  substantial  ionization  to  ensure  smooth 
plasma  ignition. 

5 .  Finite  Magnetic  Reynolds  Number  effects 

Beyond  the  zero’ th  order  approximation 
in  E,  analytical  solution  becomes  much  more 
difficult  in  the  general  case,  although  some 
additional  results  con  be  obtained  for  constant- 
area  channels.  Numerical  methods  must  then  he 
used  *o  study  the  effects  of  finite  magnetic 
diffusion.  We  will  discuss  separately  'he 
const ont-orea  case,  which  is  easier  to  analyze, 
but  somewhat  singular  in  its  properties.  Area 
variation  effects  will  be  then  studied. 

*5 .  I  The  Constant-Afea  Channel  .  from  the  work 
of  'fartinache'rand  King  it  is  well  known  that 
when  a=l  the  quas l-one-dimens lonal  conservation 
equations  admit  a  full  set  of  integrals  bv  means 
of  winch  all  variables  can  be  expressed  as 
functions  of  the  magnetic  field,  which  can  finally 
be  related  to  distance  through  Ohm's  law.  In  our 
notation,  the  first  integrals  are  Eqs.  (13)  and 
(Id’  ,  plus  ,  from  Eq.  '-4;, 


The  integration  of  Ohm’s  law  (Eq.  (16)  and 
(17))  for  the  whole  channel  then  gives 

g  -  (65) 

nj  0  ^-ub 

Eor  a  constant-  gas  model  the  conservation 
laws  can  be  combined  to  yield  a  quadratic  equation 

for  u(b):  2  (  , 

'  -(!+>:  +£  -b")G-*~-A  lb  +F. ( 1-b ) !  -  0 

2  i  -  o  o  f  to  (s^) 

The  two  branches  of  the  solution  to  (56)  are 
►he  subsonic  and  supersonic  velocities  ror  a  given 
b.  Smooth  sonic  passage  can  be  assured  by  'ore mg 
the  discriminant  of  Kq.  (56)  to  have  a  minimum 
value  of  zero  at  some  b  =  b^.  This  give,  two 
equations  relating  the  parameters  b9  and  u0*p0  t° 


►he  field  it.  Therefore,  one  can  pick  a  valu*  of 
K,  cnlculnte  u0*-p0,  solve  '56)  for  u  for  eech  b, 
and  calculate  t(g)  by  performing  the  integration 

in  (65). 

The  learnt  R  poaaible  clearly  corresponds  to 
the  maximum  of  tha  quantity, u  b,  since  otherwise 
Eq.  (65)  is  singular.  For  R  juat  slightly  above 
►bis  maximum,  the  magnetic  Reynolds  number, N»  »  1/6 
is  vei  y  large,  and  l-e  value.  -increases  ns  E  is 
increased.  An  nR*epla  for  ►.  vary  near  ►  j  n  lw 
shown  in  Fig.  I,  where  one  'an  s»*e  ’be  strong 
inlet  and  exit  current  concentration  and  flow 
m  r:e  1  era  f  Ion  reglona,  the  amn  1 1  !_evel  of  p  and  the 
predominance  of  the  back  e.m.f.  u  b  over  the  ohmic 
foil  F  u  b  over  the  body  of  the  channel.  Notice 
also  that,  despite  tbe  smallness  of  p,  its 
gradient  is  strongly  positive  n»ar  the  exit,  due 
'o  tbe  olunic  heating  in  that  region  with 
negligible  magnetic  acceleration  (small  -  bdb/df). 
Ibis  pressure  gradient  is  seen  to  decelerate  the 
flow  near  the  exit.  rhls  ‘nature  becomes  more 
gevere  as  the  magnetic  Reynolds  number  is  reduced, 
eventually  leading  to  thermal  choking,  and  is  a 
peculiarity  of  the  cons  fan t -area  geometry, 
aggravated  by  t he  somewhat  artiliciol  ^  -o 
eond  1 1 ion . 

For  an  Em  of  about  5,  Fig.  2  shows 
substantial,  deceleration  near  the  exit,  with  tho 
Vtacb  number  'dotted  ime)  approaching  unity. 

Further  reduction  to  Km  *  3.4  (Fig.  .0  leads  to 
exit  thermal  choking,  which  forces  a  mid- channel 
shock  (mathematically,  a  transition  from  'he  upper 
to  the  lower  branch  of  Eq.  (561)  and  a  subsonic 
downstream  section  ending  at  M-l.  In  this  section 
tbe  pressure  (F  is  no  longer  small.  11  is 
noteworthy  th  it  positioning  the  shock  further 
upstream  leads  to  s  lower  Em,  but,  if  h  is  held 
constant,  thw*  is  no  change  in  K  or  the  exit 
conditions,  as  long  ns  'he  flow  has  a  supersonic 
section.  A  more  realizable  condition  would  be  a 
constant  R  dimensional',  which  can  be  seen  to 
imply  b  x°Rm  -  const.,  but  numerical 
cn!rula??ons  show  only  a  secondary  effect  of  fl^ 
over  the  range  of  physical  interest. 

Eventually,  as  Em  decreases,  the  shock 
reaches  the  sonic  transition  point,  and  the  whole 
chnnnel  becomes  subsonic.  Beyond  thig  condition, 
further  reductions  of  ’bn  do  increase  E 
•  mot hrmaf ica 1 1 v ,  once  tbe  freedom  offered  by  shock 
location  is  lost,  Em  can  be  redured  in  Eq.  (551 
on.lv  bv  increases  of  c  ' .  since  bv  now  u  b<<  E 
(most 'y  ohmic  voltage),  the  ex'reme 
(electrothermal)  limit  of  Eq.  (55)  is 


The  terminal  characteristics  of  a  senes  of 

channels  with  constant  m'et  entha'py  (roughly 
i.-o,  respond  1  ag  h  -10'jK  Ai  cm!  ere  shown  in  ‘  ig.  4. 
The  leveling  of  ft,  ^  end  F  in  tbe  mixed-flow 
regime  is  clearlv  visible,  as  is  the  eventual 
Ohmic  electri^a'  behavior  (►  -vl/Eia  implies  Vw() 
and  electrothermal  thrust  variation  it’-l/Hm 
'mp'.ies  -.m2  •  Notice  'bat  ‘or  this  simple 
geometry  .  1  would  be  *be  pure  electromagnetic 
■exult,  so  that  'he  excess  s»en  in  F  at  low  Rt»  is 
tha  electrothermal  contribution  which,  as  noted 
in  Sections  4.1,  4.2,  is  non-zero  even  at  E»-><®'. 
Notice,  finally  the  quadratic  behavior  of  Er'  as  e 
-'0;  this  can  be  analytically  proven  by  expanding 
R-'ub  about  i  ta  minimum  value,  which  contribute* 


C'.CI 


the  Boet  to  Eq.  (55).  The  Rb-*  •  approximation-  is 
seen  to  provide  MO*  accuracy  for  RB  «>'Q.  Once 
«K*ln,  such  of  the  inaccuracy  is  traceable  to  the 
strong  localized  heating  peculiar  to  the  constant- 
area  geometry,  as  will  becosw  sore  clear  in  the 
foi lowing. 

5.2  Convergent-Divergent  Channel  For  a  *  a( £ ) 
the  *oaentun  equation  doea  not  integrate 
iffwued late! y,  and  the  set  of  equntiona  ll3)-(18) 
must  be  integrated  numerically  in  0°.^  1.  ror 
this  purpose,  the-  der i vat i vea  must  be  explicitly 
solved  for,  as  was  done  in  Eq.  (JJ)  for  In 

terms  of du  , 
urp 

-di  -i  '4  *  ^  (E-Gbj  <58) 

of,  a  Uq  C 


and,  of  course, 


f-(S-ub) 


Alternatively,  equation  ( *5R )  can  be  replaced 
by  algebraic  solutionof  the  total  energy  equation 
(Eq.  ',19))  for  p  -  p(u,  b). 

The  only  known  initial  condition  (at  i,-o)  is 
bfo)*l.  The  initial  velocity  u'  o'  must  be 
guessed,  and  then  p(O)  can  he  calculated  from 
pu**l  and  the  given  h  We  also  need  to  goes? 


the  parameter  fi  = 


rhe  correct  values  of 


up  and  PQ  must  Sat  is fy°( i terat ivelv )  the  two 
downstream  conditions 

'«)  Smooth  passage  through  the  ordinary  speed 
of  sound 

<b)  b( 1 )  r  C 

Notice  that  the  last  condition  replaces  the 
condition  of  magnet oson ic  throat  choking  which  was 
used  in  the  0th  order  solution.  In  ‘act,  it  is 
only  for  €  -*Q  that  strict  magnetoaonic  choking 
occurs,  Just  as  it  is  only  for  purely  Invlscid 
flow  that  an  ordinary  nozzle  flow  becomes  sonic  at 
a  throat.  In  the  presence  of  finite  magnet \c 
diffusion  effect#  (finite  Rm) ,  the  effect  of  the 
downstream  condition  b(l)=Q  can  be  felt,  althou  h 
weakly,  throughout  the  channel. 


(«)  No  current  concentration  is  seen  near  the 
inlet  or  exit,  and  the  currant  dans l t y  is  about 
uni  fora  in  tha  channel. 

( b '  Thera  is  no  positive  pressure  gradient 
(In  fact  there  is  a  negative  one)  near  the  exit, 
so  the  flow  is  never  decelerated. 

(o'  The  Mech-Alfven  numoer  (MA)  la  slightly 
less  than  unity  at  he  'hroat,  and  crosses  to  MA>1 
in  the  divergent  part  of  the  nozzle. 

These  features  generally  eonflra  the  ana’ysis 
mada  for  Riy->  despite  the  fairly  low  o*.  chosen. 
In  contrnnt. ,  the  constant-area  channel  of  fig.  1, 
for  a  comparable  Pa,  snowed  strong  e lec trot  berms  1 
features.  This  difference  clesr'.v  illustrate,  'he 
potential  of  area  contouring  for  improving  MVD 
thruster  efficiency  and,  bv  reducing  current 
concentrations,  extend  electrode  Life. 

Fig.  *5  shows  how  the  inlet  field  P(o),  inlet 
density  3(0’  and  Mnch-Alfven  number  at  the  throat 
ma  tbr. i '  evolve  as  the  magnetic  Reynolds  number 
depresses  from  infinity  to  about  3.  Clearly,  the 
salient  featuie*  of  *  be  Rmr..  limit  are 
recognizable  down  to  3,  with  only  a  slight 

increase  of  f  throughout  this  range  ithla,  of 
eourse,  implies  a  nearly  cubic  V(l) 
characteristic;.  The  throat  speed  Is  confirmed  to 
be  quite  close  to  magnetosonic  in  this  range  of 
Pm.  The  large  values  of  ('(o)  reflect  the 
smallness  of  u(o);  since  u  Increases  rapidly  away 
from  rro,  there  is  in  all  cases  a  very  repid 
density  decline,  which  will  be  very  difficult  to 
resolve  by  numerical  grids  in  multidimensional 
computations. 

6.  Comparison  to  Experimental  Data.  A  recent 
paper  by  R.C.  Jabn  and  A.J.  Kelly  provides  a  set 
of  thrust  and  voltage  data  from  a  small  M,,9 
thruster  with  no  area  dlatribution  which  is  flared 
at  first,  and  -hen  continues  for  the  greater  per* 
of  its  length  wtth  constant  cross-section.  These 
data  are  doubly  useful  in  that  they  are  presented 
in  non-dimensional  terms,  using  parameters  which 
are,  however,  different  from  ours.  The  perameters 
e-e  chosen  to  emphasize  the  importance  of  the 
Aifven  critical  ionization  phenomenon.  Thus ,  the 
choice  reference  voltage  is 


Forward  integration  of  these  equations  la 
made  difficult  by  the  existence  of  the  critical 
sonic  point  and  the  near-cn*  i  -al  magnetosonic 
point,  '’’he  procedure  followed  was  to  select 
values  of  uQ  and  and  integrate  using  the  4t.h 

order  Runge-Kutta  method  to  the  sonic  point,  where 
rapid  divergence  was  usually  encountered.  The 
value  of  F?0  was  then  modified  and  the  calculation 
re-started,  until  a  sufficiently  smooth  sonic 
passage  was  achteved.  Once  this  happened,  a 
second  divergence  was  typically  encountered  near, 
but  somewhat  beyond  the  throat,  this  forced 
repetition  of  the  whole  process  with  a  new  u0  , 
until  eventually  the  b(lj=0  condition  was  met. 


v  «  (2  b  A  v  ■ 

ref  sc 


b  -  la  ;  V  •  (•- — ' 

a "  r  a  c  u . 

c  1 


where  r  and  r  are  the  anode  end  cathode  radii,  F, 

a  l 

is  the  Tom  ration  **nf*rsrv,  f\n<4  m ,  the  ion  nan  end 
V  \«  the  Alfv«“n  critical  velocity.  The 

:i  r* 


reference  current  ,  i i 


For  Illustration,  a  convergent-divergent  area 
distribution  was  generated  with  a  -2,  suen  as  to 
obtain  a  linear  pressure  profile  between  Mach 
numbers  of  M  =0.297  and  M  =h  in  cold  Argon  flow, 
This  same  area  was  then  uaed  for  MJ’"  calculations 
wither  0,  9.1  and  0.3,  using  an  .d»al  gas  model 
with  1=  5/  3  and  h  -  0.U0328.  Th«  results  lor  c 
0.3  'Rb=3.333)  arV5  shown  in  Fig.  I,  where  several 
features  can  ba  noticed: 


We  have  converted  *he  results  of  Fig.  4  of 
this  paper  (constant  area,  Y.  '.</ 3'  to  theae  new 
parameters,  and  super  imposed  them  on  *  he  curve  of 
dimensionless  voltage  vs.  dimensionless  current  of 
uof.  0,  This  is  shown  in  Fig.  7,  which  also  con¬ 
tains  the  predictions  of  the  Aifven  critical  speed 
model.  It  cen  be  seen  that  the  fit  is  verv  good 
at  high  current  (the  un-choked  range),  but  the 
theory  under-predicts  the  linear  part  of  the 
voltage  (the  subsonic  range).  This  is  possibly 


i 


due  to  our  neglect  of  any  finite,  ionization 
effects  in  the  results  of  Ft*.  4.  These  calcula¬ 
tions  should  be  repeated  using  a  more  realistic 
(equilibrium  ionization)  gas  model.  The  important 
point,  however,  la  that  the  linear-to-eub ic  trans¬ 
ition  appears  naturally  with  no  need  to  invoke 
critical  ionization  events,  and  la  therefore  seen 
to  he  a  reflection  of  the  basic  gas-dynamics  of 
the  device. 

We  have  also  compared  the  thruat  vi.  current 
predictions  of  Fig.  4  to  the  data  reported  m 
Ref.  *J.  These  data  are  norma' i/ed  to  emphasize 
►  he  quadiatlc  dependence  on  current  at  high  I  V n 
values.  The  results  ore  shown  in  Fig.  8.  The 
theory  overpredicts  thrust  in  the  e lect romagnet  ic- 
dommated  range  'probably  due  to  some  misrepre¬ 
sentation  of  the  geometrical  effects  plus  r ne 
re gleet  of  friction,  hut  it  reproduces  very  well 
the  ►ransition  to  the  e  lect  rot  herse  t-dooi  not  **d 
range . 


7 .  Co nc fusion a  n nd  Piscu s s ion  Despite  our 
limitation  fo  quos  i  one-d  imees  iona  ‘  flews,  several 
important  features  of  seif-field  accelerated  ♦lows 
have  been  identified: 

For  reasonably  large  Magnetic  Reynolds 
numbers,  geometrical  throats  are  nearlv 
magnetoson ically  choked. 

.  Under  the  same  conditions,  the  channel 
transverse  voltage  is  dominated  by  the  local  back 
emf,  except  near  inlet  and  exit. 

Theae  two  properties  together  essentially 
determine  the  voltage  for  a  given  current  and  flow 
rate. 


Ordinary  sonic  passage  occurs  well 
upstream  of  the  throat,  and  ts  role  is  mainly  to 
determine  conditions  at  the  \njector  plate, 

Except  for  constant  area  channels,  where 
the  strong  inlet  and  outlet  current  concentrations 
produce  important  thermal  effects,  the  role  of 
ordinary  pressure  forces  in  determining  r,.ow  and 
terminal  properties  is  quite  limited. 

Approximate  factors  can  be  derived  to 
separately  quantify  the  performance  gams  due  to 
inlet  convergence  and  nozzle  divergence. 

Comparison  to  experimental  data  on 
terminal  characteristics  confirms  several  of  the 
*  rends  shown  by  this  theory.  In  particular,  the 
e  ’.act  rot  hermal  -electromagnetic  transition  is  well 
predicted,  with  no  need  to  invoke  aulation  or 
critical  ionization  effects. 

Many  of  these  properties  are  likely  to  be  of 
■se  In  mul  t  id  1  mens  tone  i  Bicx'e  tng  of  wp»j  flow*. 
v,,r  instance,  approximate  f’uw  fields  can  be 
'••imputed  much  wore  ensilv  by  neglecting  axial 
pressure  forces  wnile  retaining  them  in  the 
"rnnsverae  directions  to  balance  pinching  magnetic 
fo»cM;  difficulties  associated  with  ‘ransonic 
f'«)w  are  thus  obv'a*ed.  It  \s  also  usetul  to 
recognize  the  dominance  of  the  bark  emf  over  ohmic 

provided  current  concentration  layers  are 
treated  separately.  Finally,  such  current 
'  onrent rat  ions  nan  be  minimized  by  the  same  sort 
of  careful  aerodynamic  shaping  as  exerc.sed  by 
ordinary  nozzle  designers. 


To  conclude,  it  must  be  recalled  that 
throughout  this  analysis  we  have  assumed  no  net 
axial  current,  and  hence  no  gross  transverse 
nonuni for-mity.  While  this  is  a  realizable 
condition,  many  experimental  accelerators  are  in 
fact  strongly  affected  by  axial  current  and 
pinching  magnetic  forces,  so  that  our  results  have 
only  a  qualitative  validity  in  such  cases;  in 
particular,  the  dramatic  electrical  effects  of 
anode  depletion  ( ’'onset'*)  are  not  predictable  by 
this  type  of  tneory. 
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Abstract 

A  two-dimensional  numerical  model  has  been  developed 
in  order  to  analyse  electro- magnetic  p'asma  accelerator* 
also  called  Self-Field  Magneto- Plasm  a- Dynamic  Thrustere. 
This  model  uee*  a  Magneto- Hydro-Dynamic  description  of 
the  gas  considered  aa  a  'ul!y  ionised,  isothermal  peasma, 
and  take*  into  account  the  Hall  effect  (non  linear  conduc¬ 
tivity)  and  the  interaction  between  the  magnetic  field  and 
the  fluid  dynamic*  of  the  plaema.  The  eyetem  of  equation* 
i*  diecretiied  into  finite  volume*,  and  i*  *o!ved  by  a  Newton- 
P-aph»on  scheme.  Reeult*  from  the  MHD  model  were  calcu¬ 
lated  for  a  maee  flow  rate  of  6  g/«  o'  Argon  and  'or  current* 
up  to  ten  kilo-Amperes. 


Introduction 

During  the  paet  year*  the  development  o'  high  current 
plaema  accelerators  (also  called  Magneto- P'.aema- Dynamic 
"Truster*  )  has  been  aimed,  among  other  thing*,  at  im¬ 
proving  the  efficiency  of  energy  con/ersion  from  electrical 
to  kinetic  energy.  Several  researcher*  have  observed,  that 
the  efficiency  increase*  with  the  current,  but  at  som-  va  ns 
a  transition  occur*  to  an  unstable  and  deetructive  regime  of 
the  discharge  (called  onset  of  instability  in  the  literature), 
when  tni  ratio  ■—  exceed*  a  critical  value  which  depend* 
on  the  nature  of  the  propellant  and  the  thruster  geome¬ 
try.  Some  researcher*  3’  'U'  have  proposed,  a  theory  for 
the  transition  to  an  unstable  behavior  based  on  the  deple¬ 
tion  of  the  anode  from  current-carrying  electrons,  the  for¬ 
mation  of  a  starvation  layer  along  the  anode  characterised 
by  current  den*itie*  in  exceM  of  the  random  thermal  flux  o' 
electron*,  and  the  reversing  of  the  potential  drop  between 
the  the  plaema  and  the  electrode.  In  order  to  study  the 
'"ormat.on  of  thi*  *tarvation  layer  one  needs  to  take  the  geo¬ 
metrical  properties  of  the  discharge  into  account,  or  in  other 
word*,  to  look  at  a  two-dimensional  mode!  of  the  discharge. 
Eventually  thi*  model  can  be  used  to  find  opt’mal  shape* 
that  maximi**  certain  parameter*  like  efficiency  or  electrode 
life. 
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Various  attempt*  have  been  made  to  develop  numerical 
two-dimensional  modei*  of  MPD  discharge*  u*ing  finite  el¬ 
ement  or  finite  difference  description*.  1.  fCmura  et  ai.  17' 
have  presented  a  two  dimensional  modei  for  the  electromag¬ 
netic  field  (self  field),  including  the  Hall  efiect.  The  model 
was  primarily  designed  to  calculate  the  current  distribution 
near  the  cathode.  They  later  extended  their  mode!  27'  to 
systems  that  work  with  an  external  focussing  magnetic  field. 
These  model*,  however,  were  not  de»:gned  for  an  accurate 
modelling  of  realistic  thruster  geometries,  and  did  not  in¬ 
clude  the  effect  of  the  convection  of  the  electromagnetic  field 
by  the  plasma.  This  last  limitation  was  partially  removed 
in  a  later  paper  (28:  where  a  simplified  calculation  of  the 
flow  was  made  and  coupled  with  the  magnetic  field.  Ao  and 
Ti; iwara  (1  have  attempted  to  develop  a  one-fluid  Magneto- 
Hydro-Dynamic  mode’  of  the  plasma  where  the  equations 
d.esc-lb'ng  the  motion  of  the  plasma,  considered  as  a  neu¬ 
tral  conductive  fluid,  are  coupled,  with  the  magnetic  fie.d 
through  the  Lorents  force  and  the  ohmic  heating.  However 
this  mode!  was  Limited  to  a  very  simple  geometry. 

Tt  was  decided  to  develop  a  two-dimensional  numerical 
ni"d*  based  on  a  MHD  description  o'  the  plasma  in  the 
collision  dominated  reg'ine,  ,n  order  to  ca.culats  th*  steadv 
state  d  stribution  of  the  electro-magnetic  field  and  the  mo¬ 
tion  of  the  plasma  ;n  the  accelerator.  We  focused  our  atten¬ 
tion  on  the  phenomena  observed  at  the  electrod.ee,  especially 
at  the  anode,  :n  an  efiort  to  detect  the  onset  o'  the  turbulent 


The  system  o'  partial  differentia!  equation*  is  discretised 
using  a  finite  volume  technique  used  extensively  in  Compu¬ 
tational  Fluid  Dynamic  algorithms.  The  system  of  equa¬ 
tion*  was  solved  using  a  Newton- Raphson  method,  derived 
from  an  initial  calcu.ation  o'  the  magnetic  field  alone  for 
which  this  method  works  well  due  to  the  dominant  elliptic 
behavior  of  the  equations.  The  method  o'  discretisation  of 
the  equation*  (finite  volume  method.'  as  well  as  th*  idea 
of  the  addition  of  an  artificial  vi»co«ity.  was  derived  from 
Jameson  and  Yoon  15’.  A  confirmation  of  the  possibility 
o'  the  Newton- Raphson  method  was  found  in  the  lecture  o' 
Jespersen  ’Hi’,  as  w»l!  as  in  th#  doctoral  these*  of  M.  Giles 
fi;  and  M.  DreU  ’7‘ 
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Physical  Model 


The  equation*  describing  the  model  are  derived  from  the 
general  Magneto-Hydro- Dynamic  approximation  which  in¬ 
clude*  the  following  equation*:  conservation  of  dim,  mo¬ 
mentum  and  energy  in  the  inviacid  fluid;  the  Ohm’i  '.aw 
(derived  'rom  the  electron  momentum  conservation);  and 
the  quasi-static  Maxwell's  equations. 

dip  -t-  V(/ra)  =  0 


9,pn  4-  Div(puu  —  S  -  M)  =  0 

Where  5  is  the  pressure  stress  tensor,  and  .Vf  Is  the  Maxwell 
stress  tensor. 

9,(Em*,  +  Ei  +  EK)  +  V  ({F,  +FK  +  P'!n+T— ? +a)  -  0 

P’0 


Where  East,  Ei,  Ek  are  respectively  the  electromagnetic 
energy  (  B1/2po),  the  internal  energy,  and  the  kinetic  en¬ 
ergy,  per  unit  volume;  P  is  the  pressure;  *  >*  the  Poynt- 


ing  vector;  q  is  the  heat  conduction.  The  snc-gy-equation  is 
usually  completed,  bv  a  state  equation  relating  the  pressure 
to  the  internal  energy,  for  Instance  ; 


P  =  (y  -  1  )E, 

The  generalised  Ohm’s  law  can  be  written  ae: 

E*«xfl=-J  +  -l-  J  x  B  -  V-- 

<7  cn,  m. 

The  Maxwell's  equations  reduce  to: 


V  x  E  =  -9,  3  VB  =  0 

V  x  B  =  MqJ 


Ws  nave  assumed  that  the  fleld  is  quasi-static  and  that  ths 
p  aema  .*  approximately  neutral. 

Let  ue  now  look  at  the  possible  simplifications: 

Geometry:  The  flow  and  the  electromagnetic  fleld  are 
ax.8ymetric.  There  is  no  externally  imposed  magnetic  fleld. 
The  Magnetic  field  .s  therefore  asimuthal. 


The  Ohm’s  Law  and  the  electromagnetic  equations  are 
combined  in  order  to  eliminate  E  and  J.  One  obtains  an 
equation  in  B  that  describee  the  difiue  on  ar.d  convection  of 
the  magnetic  field  in  a  conducting  medium: 


_  ,  „  V  x  B  fV  x  B'  x  B  VP, 

fMj  en,  O', 

Rewriting  the  equation  in  non-dimensional  form. 


-d.-B*  4-  V*  x  (o’  X  B")  = 


x  (  V*  *  B* 


'V*  X  bo  X  B‘ 


2  n* 


introduces  three  non-dimensional  numbers:  Si  =  jf—  for 
the  magnetic  diflueion  due  to  conductivity;  Si  —  for 
the  non  linear  Hall  e^ect;  S»  =  for  ths  electron  dif¬ 

fusion.  Pm  =  pqCoMoLo  s  the  magnetic  Reynolds  number, 
based  on  the  d.stance  between  ths  electrodes;  Ht  .»  the  Hall 
parameter; 


Bo7q 

- =  U/..T 

rn.o 


where  w„  is  ths  electron  cyclotron  angular  frequency  and 
r,  =  1/m,  is  the  electron  collision  time  .  0  is  the  plasma-;?; 
ratio  of  ths  pressure  to  the  magnetic  pressure, 


P  = 


2/ieP. 


The  values  of  thee#  non-dimensional  numbers  give  an  indica¬ 
tion  of  ths  relative  importance  of  each  phenomenon:  Inside 
the  channel  ,V(  63  0.5,  Si  rj  0.0,  S,  cs  0.03.  Therefore  we 
take  the  Hall  effect  and  the  convection  into  account,  but  we 
neglect  the  effect  of  the  electron  pressure  gradient  on  the 
conductivity.  (Notice  that  a  magnetic  Reynolds  Number 
based  on  the  channel  length  would  be  of  order  5  to  10). 

Constant  temperature  approximation:  T.n  a  plasma  like 
the  one  considered  here,  the  energy  equation  is  the  most 
difficult  to  justify.  At  high  temperature  (around  12500 
K)  the  gas  is  strongly  ionised  and  radiates  energy.  These 
processes  are  not  necessarily  in  equilibrium  and  therefore 
cannot  be  accurately  described  by  the  energy  conservation 
equation  used  In  the  fluid  theory  at  ’ow  temperatures.  As  in 
meet  V_HT  modex  ths  energy  conservation  equation  w:li  be 
dropped  and  rep  \ced  w;th  the  assumption  that  the  plasma 
:*  roughly  isothermal.  This  assumption  can  be  justified  by 
the  eBect  of  'on  « at  on  wh  ch  tencs  to  oppose  any  change 
n  temperature  '•»  a  d.sp.arement  of  the  on  tat. on  equi':  u- 
•  'j'j,  rr  the  thermal  conductvty  of  the  plasma,  and  by 
fe  value  o'  the  pressure  relative  to  the  raagnet.c  pressure 
’’’h#  non-d . mene’on ,v  number  deecrib.ng  ths  ba.enc* 
between  kinetic  energy  and  heat  conduction  (equivalent  of 
the  Reynold*  number  'or  the  energy  equation)  is: 


n  Up0r;r 

c,  =  —  - — 

k1  1 

Its  value  goes  from  «j  200  in  the  cathode  jet,  to  tf  1  along 
the  external  •ur'ace  of  the  anode. 


Conductivity:  In  the  Cpitier-Fvrm  aproximation  the  con¬ 
ducts  ty  o'  a  'ul'v  on  r"d  plasma  .s  proportional  to  T 1 
Consequently  w#  have  a  constant  value  'or  the  electrical 
conduct. v  tv  <7  =  ?a,XSi/m. 

Even  though  the  enrgy  balance  is  not  written,  the  'act 
that  one  part  of  the  »  »c*.r-cal  work  E  J  s  dies  pated  ohmi- 

(  J 

cally  as  ■-  s  ;mpllc:t  in  the  rest  o'  the  formulation,  and 
therefore  1  meaning'ul  thruster  efficiency  can  b*  computed. 
What  is  igno-ed  is  the  ultimate  'ate  of  ’he  ohmic  h»at  ar.d 
ts  apport  nnment  to  onisat.on,  heat  conduction,  rad. at. on, 
etc. 


2 


Since  ionliation  processes  are  not  covered  in  the  model, 
the  p'.umt  n  considered  fully  ionised  when  it  enter*  the 
inlet  of  the  accelerator. 


Although  the  real  Reynold*  number  !•  of  the  order  of  ’00 
to  200,  the  plaema  le  coneidered  inviecid,  which  mean*  that 
second  order  derivativee  appear  only  in  the  electromagnetic 
equation. 


The  final  set  of  equation*  after  simplification#  1*  then: 

1  , 


9,(pu.)  j9,(rpu,']  =  0 


J93  l 

9, {pH,  u,  +  P  +  - — )  +  ~9,(rpu.  ur)  =  0 
-IM>  r  ' 


9,{p%.H ,)  +  -9,(r[pu,U,  P  +  - — " 

r  2  in, 


B'  -P+-?- 

-_M  +  - fjLL  =  0 


P  =  pRTa 
Ed 1  =  0 


/ 


E.  =  —  I  8,rB - - — 9,B  —  u,B 

IM)<T  r 

E,  =  -  —  8.B - 4-  u.B 

Ho<r  Morn,  r 


The  boundary  condition*  were  taxer  a*  follow*: 

The  tangential  electric  field  i*  »ero  on  the  electrode*.  The 
electrode*  include  the  conducting  surface*  of  the  anode  and 
the  cathode.  The  voltage  drop  aero**  the  p'.asmvsheath, 
which  extendi  over  a  few  Oebve  length*  inside  the  plaema, 
w;:_l  be  neglected.  Coneequente’.y  the  voltage  drop  acroee  the 
channel  will  be  undereet  mated  w.nen  the  anode  1*  etarved, 
at  high  current  deneitiee. 

The  entrance  surface  (the  back  plate)  i»  an  insu'.ating 
surface  which  carriee  the  injector  holee.  It  1*  assumed  that 
the  gae  n  injected  uniformly  acroea  the  back  plate,  which 
would  be  the  case  with  a  ooroue  back  plate  or  with  a  olate 
conta,n:ng  a  large  number  of  injector  hole*.  A  transition 
between  the  gae  *tate,  charactcmed  by  a  low  temperature 
ard  a  low  conductivity,  and  the  p'aema  state,  where  the 
gae  ie  ionised  and  ie  a  good  conductor  of  electricity,  arisee 
or.ee  the  gae  ie  in  the  annular  chamber.  Thie  :onisation 
takes  place  in  a  region  eioee  to  the  back  plate  surface.  The 
cactilation  starts  after  the  Ionisation  'aver  assuming  that 
the  gae  hae  reached  a  uniform  t*mp*ratur*  of  UlOO  K,  at 
whicn  the  elect  ncal  conduction  ir  «u*ficent'.y  Vgh  to  consider 
the  gae  a*  a  plaema.  Acros#  the  entrance  plane  the  magnetic 
fie.d  B  —  Ho / / 2 v e .  The  dynamic  variable*  p,pn,,p u, 
are  specified  a.ong  the  entrance  surface  :  p*i,  =  0;  pH,  is 
specified  by  the  mas*  flow  rate: 


pH,  = 


d,nl«ti 


The  existence  of  a  steep  density  gradient  immediately  ad¬ 
jacent  to  the  back  plate  make*  it  very  d:*fi.cult  to  reeolve  this 
layer  numerically.  Fortunately  one  of  us  has  shown  22  that 
the  detain  of  the  Inlet  gTviient#  have  only  a  »ma!l  influence 
downstream.  In  our  calculation  we  assign  an  arbitrary  value 
to  p. 

On  the  outer  boundary  of  the  computational  domain  (i.s. 
far  into  the  piume)  B  «■  0.  Thi*  condition  means  that  there 
I*  no  current  going  outside  the  boundary  of  the  calculation, 
which  is  therefor*  an  insulating  surface.  There  are  no  gra¬ 
dient*  of  p,  pu,,  p%,  along  the  direction  parallel  to  the  flow. 

On  the  axis  of  symetry  B  —  0. 


Numerical  Method. 


The  system  of  equation*  i*  solved  by  the  Newton- 
Raphston  method  using  the  technique  of  Finite  Volumes 
for  the  discretisation  of  the  equations  of  motion,  and  a  tech¬ 
nique  similar  to  the  method  of  finite  volume*  for  the  mag¬ 
netic  field.  Artificial  viscosity  is  added  to  the  equation*  of 
motion  in  order  to  avoid  non-linear  instabilities  at  places 
of  strong  density  gradients.  The  conservation  equations  are 
discretised  according  to  the  method  of  finite  volumes:  The 
computational  space  is  div'ded  into  20  x  41  quadrilaterals. 
-T’he  gr id  Is  shown  in  figure  1.  The  variables  (p,  pu,,  p%,,  B) 
are  aoproximated  by  their  space-average  on  each  quadrilat¬ 
eral.  ''.'he  equal  "ns  are  written  In  the'-  Integra,  form: 


where  the  the  tensor  pn u  —  E  —  M  ha*  the  following  physical 
component*  in  cylindrical  coordinate*  [r,r,  ■. 

r  ,  h 

PH.H.  -t-  P  PH,  U,  0 

_  n ) 

pH,  U,  PH, Hr  P  +  y—  0 

0  0  P  -  5®-  | 

Since  the  problem  is  axisymetrical,  the  volume  elements 
are  obtained  by  rotating  the  auadr.lateral*  that  define  the 
gnd  cells  by  an  infinitesimal  angle  Ai t>  around  the  axis  of 
symetry.  ’’’he  sur'ace  integral  around  each  element  :s  bro¬ 
ken  into  six  integrals,  on*  for  each  'acet.  Kach  o'  them 
is  discretised  bv  assuming  that  th*  fiuars  ar*  constant  on 
each  facet  o'  th#  boundary,  ’'or  '‘■s  electro-magnetic  field 
one  uses  a  d  fierent  integral  (line  'ntegra.  over  a  cioeed  con¬ 
tour): 

/  Edl  =  0 


'rhe  contour*  are  defined  bv  th#  same  quadrilateral#  as 
above.  vor  each  iteration  the  Newton- Raphson  algorithm 


n 


invert*  *  matrix  of  order  3200  which  represent*  the  *y*tem 
of  partial  differential  equation*  linearized  at  the  current  es¬ 
timate  of  the  solution  vector.  Thi*  take*  about  10  minute* 
on  a  microVax.  On#  need*  about  10  iterations  to  attain 
convergence. 


Experimental  Background 

Experimental  studies  have  identified  three  modes  of  op¬ 
eration  for  the  plasma  accelerator:  electrothermal  mode, 
for  which  the  dominant  acceleration  mechanism  is  the  ex¬ 
pansion  of  the  heated  ga*  acroee  the  thermodynamic  pres¬ 
sure  gradient;  electromagnetic  mode,  for  which  the  d.om- 
nant  acceleration  mechaniem  la  the  expansion  of  the  gas 
across  the  magnetic  preseure  gradient;  finally  the  unstable 
regime  which  set  on  when  the  critical  ratio  L-  is  exceeded. 
For  argon  a*  propellant,  Huge!  111'  gives  an  empirical  value 
for  the  limiting  current, 

—  «  2.5  1010  (A1  kg-1  e)  , 

fr\ 

which  givee  a  current  of  the  order  of  12  VA  for  rn  =  fl  g/s. 
More  recent  work  at  Princeton  and  elsewhere  ha*  raised 
thi*  by  factor*  of  up  to  3,  depending  on  the  geometry  and 
the  injector  configuration.  An  oneet  prediction  based  on 
tns  total  starvation  of  the  anode  was  proposed  by  Baksht, 
Moiihe*  and  Rybakov  [3], 

P  m3  _  945  ekT  LH 8 
r h*  8  <Tnf,  d * 

which  lead*  to  a  critical  value  of  19  kA  for  our  geometry. 

An  explanation  which  does  not  include  the  Hall  effect 
ha*  been  proposed  by  Lawless  and  Subramanian  19  ,  based, 
on  the  mechanism  of  an  excessive  back  EMF  induced  by  a 
strong  convective  effect  t>  x  B  in  the  opposite  direction  of 
the  applied,  driving  field  R.  Neglecting  the  Hall  effect  they 
obtained.,  from  a  one-dlmemsional  analysis,  an  expression 
for  the  limiting  current: 


/*3  a* 

- 8.52  — 

m  tMs 


(A1  kg“  .) 


where  a*  is  the  speed  of  sound  at  the  choking  point  and  no  is 
the  permeability  of  the  vacuum.  This  givee  a  current  of  the 
order  of  20  kA  for  a*  =  10*  m/s  and  m  =  6  g/s.  However 
the  physical  basi*  of  their  argument  is  open  to  criticism  22 ' . 

Assuming  that  the  magnetic  field  is  confined  between  two 
coaxial  cylindrical  electrode#  and  does  not  extend  outside 
the  exit  plane  of  the  thrusts'-,  one  can  integrate  the  effects 
of  the  Maxwell  strees  tensor,  and  one  can  predict  the  thrust 
due  to  the  electro- magnetic  effects  by  the  formula: 

r.„  =  ^inr-* 

sr  rc 

1 1 1  Newton  at  lOkA  for  r*  =  3)  The  ratio  gives  the 
exit  velocity  c,. 


The  voltage  is  expected  to  be  linear  with  the  current  for 
low  current#  'Joule  heating',  and  to  Increase  at  a  higher 
rate  'or  high  currents  as  the  field  accelerate#  the  p!a*ma 
and  compete#  against  the  induced  back  EMF,  showing  a 
cubic  dependence  with  the  intensity. 


Results 

We  have  nun  fully  coupled  calculations  for  a  mass  flow  rate 
of  6  grams  per  second  and  for  various  electrical  currents  be¬ 
tween  9  and  19kA.  Uncoupled  calculation#  (i.e.  calculating 
the  electromagnetic  field  with  a  froxen  low  field)  have  been 
obtained  for  currents  between  0  and  100  1<A. 

A  maximum  of  the  magnetic  field  is  observed  along  the 
cathode.  Figure  2  shows  the  contoun  of  constant  rB.  For  an 
axusymetrical  geometry  these  line#  also  represent  the  lines 
of  current.  The  influence  of  the  magnetic  field  on  the  con¬ 
ductivity  tensor  (Hall  effect),  is  teen  in  the  way  the  current 
line*  concentrate  at  the  tip  of  the  anode  and  at  the  root 
of  the  cathode.  As  the  current  increases  this  effect  becomes 
more  pronounced.  This  effect  can  also  be  obeerved  on  figure 
3,  which  shows  the  current  density  at  the  electrodes. 

Figures  4  and  3  show  the  density  contour*  for  0  and  lOkA. 
For  lOkA  (figure  5)  the  density  shows  strong  gradients  and 
large  fluctuations.  At  the  inlet  the  density  drop#  abruptly 
over  a  few  cell*.  The  compression  of  the  plasma,  under  the 
influence  o'  the  Lorenti  force  (inward  pinching  effect),  can 
V-  observed  a  ong  the  cathode  and  along  the  axi»  of  symetry. 
Fgure  0  shows  an  enlargement  o'  the  density  plot  for  lOkA, 
when  one  observe#:  the  -ompression  o'  the  plasma  along  the 
cathode;  the  strong  rarefaction  that  follows  the  expansion 
of  the  plasma  around  -..he  cathode  comer,  which  reaches  a 
minimum  of  ’.5  10J'J  m_s;  the  recompression  along  the  axis 
o'  svmetry,  where  the  density  reaches  a  maximum  of  108. 
19J0  rn-5  at  a  distance  of  ’.  cm  from  the  tip  of  the  cathode. 
On  the  anode  side  the  plasma  expands  monotonical’v  from 
33.  10JO  rn - '  to  0.15  lO30  '■'igure  7  shows  the  density 

across  the  thruster  a.ong  three  different  rows  of  cel's.  It 
shows  the  considerable  increase  in  density  in  the  cathode 
jet,  downstream  of  the  cathode  along  the  axis  of  symetry. 

A  plot  of  the  velocity  vector  obtained  'or  lOkA  (figure 
8  '  shows  a  strong  acceleration  in  the  surface  layer  that 
surrounds  the  cathode.  The  existence  of  this  iayer  can  be 
n'er-ed  f'om  the  Ygh  values  o'  both  the  electric  fie'd  and 
the  current  density  along  the  cathode  (about  an  order  of 
magnitude  higher  than  a*,  the  anode). 

The  dot  of  the  Hall  number  ff-  (figure  9)  shows  con¬ 
siderable  gradients  on  the  anode  and  cathode  tips.  These 
gradients  are  related,  to  the  electron  density  gradients  cre¬ 
ated.  bv  the  expansion  of  the  plasma  as  it  follows  the  elec¬ 
trode  surfaces.  The  Had  effect  (Tentorial  conductivity)  :s 
likely  to  be  responsible  for  the  computational  inetab.lities 
encountered  at  higher  rur'ents  along  the  electrode  surfaces. 


The  thrust  was  computed  for  various  current*  and  is  lifted 
in  table  2.  The  thru*',  increase  due  to  the  current  ;*  appar¬ 
ently  very  modest  (9.8  N  at  lQkA),  but  is  comparable  to 
the  thruit  due  to  the  Maxwell  »tree*  tensor  predicted  by 
the  formula  : 


(11  N  at  lOkA  for  (*■  =  3  ).  Thu#  the  thrueter  ii  etill  oper¬ 
ating  at  this  condit.on  mainly  in  the  electrothermal  regime. 

Conclusions 

The  calculation  of  the  MHD  mode!  was  done  for  three  val¬ 
ues  of  the  current:  9,  5.VA  and  ’OkA.  The  mode!  does  not 
converge  above  10k A.  The  probable  cause  of  the  computa¬ 
tional  instability  is  the  excessive  value  of  the  Hall  parameter 
along  the  surface*  where  rarefact  on  is  present.  Above  a  cer¬ 
tain  value  of  the  current  this  rare'act:on  will  be  increased 
by  the  action  of  the  Lorents  force,  leading  to  a  failure  of  the 
model  along  the  boundaries. 

The  sam«  cause  may  be  responsible  for  the  onset  of  insta¬ 
bilities  observed  in  experimental  devices.  (The  onset  cur¬ 
rent  expected  from  Huge!’#  formula  is  12kA,  the  value  from 
the  Baksht  et  al.  is  19kA,  and  the  value  from  Lawless  and 
Subramanian  is  20kAj. 

The  artificial  damping  introduced  in  the  equation  of  mo¬ 
tion  in  order  to  obtain  a  reasonably  well-oosed  problem 
modifies  the  solution  bv  an  amount  that  reaches  8  percent 
of  the  momentum  fiux  (for  10k A'.  The  largest  contribution 
to  this  error  appears  at  the  sharp  turn  around  the  anode 
corners. 
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Symbols 


P  F'us).B<ta  'oreesure/magnet.c  pressure 


Spsciflc  heat  ratio 

A 

Heat  conductivity 

A*o 

Permittivity  of  vacuum 

a* 

Speed  of  sound  at  chocking  point  in  the 

P 

Mass  density  (  kg  m“*  ) 

theory  of  Lawless  and  Subramanian 

<7 

Conductivity  (  Si  m_1  ) 

B 

Magnetic  field  (T) 

% 

Electron  collision  time 

c« 

Exit  velocity 

<t> 

Axisymertic  coordinate 

Ct 

Non  dimensional  parameter 

W„« 

Electron  gyrofrequency  f angular 

In  the  energy  equation  (Clausius  Number) 
d  distance  between  electrodes 

t  Charge  of  the  proton 

E  Electric  field  (  V  m-1) 

Ebm  Electromagnetic  energy  per  unit  volume 
Et  Internal  energy  per  unit  volume 

Ek  Kinetic  energy  per  unit  volume 

H  Average  perimeter  of  the  Inlet 

Hall  parameter 

/*  Current  at  the  threshold 

of  the  turbulent  regime. 

J  Current  density  (  A  m~’  ) 

k  Boltsman  constant 

L  Average  length  of  the  ehanne' 

m,  Volecu’ar  mass  of  the  gas  (  kg) 

M  Maxwell  stress  tensor  (  N  m-’) 

m  Mass  How  rate  of  the  propellant  (kg  s-1' 

rn,,rnr  7.-  and  r-momentum  (=  pu,,pu,) 

n.  Density  of  electrons  '  m~*  ) 

P  Pressure  (  V  m'1  l 

q  Heat  flux  (  W  m-’  ) 

R  Gas  constant  (  J  kg-1  K~l  ) 

R-m  Magnetic  Peynolds  number 

o’  Stress  tensor  '  N  m-J) 

T  Temperature  (  K  ) 

a  Velocity  (me-1  1 

iu  Diad  (u,«,  :n  cartesian  coordinates' 
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& 


/oV  A  ' 


n 

Cathode  radiue 

A 

T. 

2.  cm 

• 

J  Anode  radiue 

r. 

6.  cm 

V  Cathode  lenj^th 

6.  cm 

|  Anode  length 

12.  cm 

£*•  Inlet  deneity 

Pml.t 

2.21  KT4  kg  m-® 

fc'1  Inlet  preeeure 

R.nl.t 

1144  Pa 

Inlet  velocity 

Uz.nl., 

2700  m/e 

M  Inlet  X-eec.  area 

A4*/#i 

10-3  m’ 

t_-  Temperature 

To 

1.25  104  K 

p.'  Conductivity 

tto 

3.8  10s  Si/m 

l 

*"  Vase  Sow  rate 

m 

*  «/« 

£  Current 

1 

(0.,  5.,  10.)  kA 

i 

y  Gae  constant 

R 

0.416  kJ/kg/K  ! 

'.0. 
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F* i r«  2:  Current  ?or  /  «*  10^/4  f'.'nVu'nt'on) 
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Figure  7:  Density  across  the  channel  'or  I  —  '.OkA  (Calculation)  'a'  6th  row  o'  ceils  lns.de  the  channel';  V  14th 
ce..s  (end  of  the  cathode);  (c)  7 J t h  row  o!  cell*  beyond  tne  t.p  o'  the  cathode'. 
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Abstract 

"Vo  annular  magnet oplaamadynaaic  ^fPD > 

thrusters  were  constructed  and  operated  in  a 
regime  where  anode  starvation  processes  .unit  the 
thrusters'  per  fo  nuance .  Magnetic  ‘‘ie_d  and 

Moat  mg  potential  measurements  were  conducted  to 
determine  the  uperat tona.  differences  between  a 
channel  with  constant  intere.ee *■  rode  separation 
and  a  channel  with  an  initial  cons  tart  v*paraf  ion 
and  a  divergent  exit.  "ha  channel  with  Me 
constant  electrode  separation  shown  evidence  o+' 
strong  concentrations  of  current  a*  boM  inlet  and 
exit,  along  with  strong  indications  of 
accompanying  electrode  erosion.  In  the  divergent 
section,  current  densities  were  gr*at  1  v  diminished 
with  aecreaaed  evidence  of  ablation.  Evidence 
tends  to  indicate  that  e.ec^roce  erosion  is 
unavoidable  for  a  starved  anode:  however,  through 
spatial  variation  of  the  intellect roce  spacing 
local  ohmic  heating  can  be  controlled  to  minimize 
ablation  of  the  otrier  components,  sueK  as  the 
cathode  and  backplate  msu-atoe. 

Noaenc lature 

a  aagnetoacoust ic  speed 

"  magnetic  induction  field 

?  electric  field 

s’  v.h’-’ust 

J  total  current 

current  density 
in  mass  flow  rate 

R  gas  constant 

r,  anode  radius 

rc  cathode  radius 

f  plasma  temperature 

a  mean  mass  ve.ocitv 

?  Hal  1  parameter 

v  specif. c  heat  ratio 
m*  permeability  of  free  space 
3  plasma  sea. nr  conductivity 

Introduction 

The  goal  of  th«*  designer  of  an  *PD  arejet  is 
to  obtain  desirable  overall  per  ''onnn«  e .  "*0 

achieve  this  goal,  1 1  is  use‘‘u.  ro  ■nanipulere 
internal  distributions  of  current  tens: tv  ind  mass 
flow.  In  a  previous  paper,  Martinez  K*mcne?  and 
Reimerdinger  derived  an  approximate 

two-dimensional  hydrodynamic  YPO  arc jer  mor.ei 
which  indicated  that  *  he  performance  and  "he 
internal  distributions  of  the  previously  described 

This  work  was  supported  by  the  Air  force  Office  of 
Scientific  Research  under  grant  number 

APJSR-8*-0119D. 


quantities  can  be  manipulated  by  an  axial 
variation  of  1  ha  inter#>rtr'iy  separation.1  Ms 
P*per  presents  prei  immarv  lata  '  r  )tn  an  experiment 
based  on  the  analysis  from  reference  . 

MpQ  arcwe‘s  can  be  distinguished  from  other 
electric  thrusters  by  their  use  o‘  *  he  Jxfl  Lorenz 
force  as  tneir  primary  source  of  thrust.  An 
approximate  overall  thrust  equation  snows  ‘'his 
future  by  the  quadratic  behavior  of  the  Mrust 
with  the  applied  current  and  a  weak  dependence  ,.,n 
geometry  1 


for  a  given  thruster  geometry  and  mass  'low  rate, 
the  efficiency  a. no  increases  w  M  ut rent  because 
of  the  transition  rroa  e  .e<  t  ror  hema  .  acceleration 
due  to  the  ohmic  heating  of  ‘he  plasma  to 
primarily  an  e iect romagret ic  acceierat ion . 

.n^ortunatel v ,  an  operational  1  .ait  called  'onset  ' 
is  reached  where  b^M  erosion  o'  the  Mruater  ind 
overall  voltage  Mve  be**n  n  o‘ed  to  mc^ease 

■'ama‘‘ca  ly.  ""he  ‘*rm  onset  originated  *>r^s  '  •  :e 
*  of  Tieghner‘7  tec^j:  nai.  v<>  N.’i* 

.  a'  -»*e  whit;1  i-*'..*  v  grow  nap  , 

tj.’tA.ph  lufnet  onset  ns  Mat  ev*l  w*are  vu  ‘ape 
< -a*.  1  '  '  at  1  r>n  r  *Sa  h  .OX  >'  ‘he  average  ►en*ir'nl 
vo.’ngn.*  fn  a  */:ven  geometry,  *  he  onset  m  •  r 
s  n.es  approximate  v  with  the  pernae'**'-  J*  m  nj;d 

•  h  -is  corre  stes  With  ‘hi-  exhaust  ve  1  v*  1 1  y  .  •  •  1 

'n  --‘Vrericw  ..  *  he  »pp**a' ance  of'  onset  is 

attributed  *0  an  *\n*'de  starvation  phenomenon 
-  a 'used  bv  axin’  de  ♦'  \  e«.  *  \  **«r  o‘*  ’he  local  'urrent 
vector  *' y  ‘he  M.*  effect.  In  figure  ,  "he 
;  r  en  t  j  firm,  wh  ’  ■  h  is  orthogonal  boM  ‘o  ‘he 

ic  .  M.«uc‘:on  e  eld  and  ‘he  current  density 

•  e.  tor,  has  s  .  moonen  t  directing  the  plasma 
‘  wards  ‘he  «  at  hoc e.  \s  *  he  plasma  adjacent  ‘ o 
‘he  anode  cecomea  r a relied,  deflection  of  the 
current  ve*  *  or  is  exacerbated  ns  the  wal. 
.>nramn'«r  g  r  ■  ■ws  due  to  'he  o«  a !  -'ecreasa  m 

sw.na  ‘v.  "he  *-egion  adjacent  to  ‘he  <\node, 
ihi.reaaes  .  n  vo.'age  "  o  at‘rac‘  electrons  ar.d 
ensure  current  fli.ag.  Eventual  lv  the  magnitude  o ♦' 
Mis  potential  grows  large  enough  for  e. metrical. 
b*"e»ikd<jwn  ‘*>  oi:cur  wh.ch  a.  lows  concentrated  arcs 

•  o  form  on  ‘he  anode.  rhis  ef*‘ect  is  intrinsic  to 

arejet  operation  un  .ess  an  external  electric 

•  :eld  is  applie.j  to  halt  the  growth  o*  axial 
current . 

In  addition  to  plasma  depletion  at  the  anode 
at  high  currents  and  low  m#iss  rl«jw  rates,  ‘ha 
arejet  can  operate  at  a  significant  magnetic 

•  e y t  1 « »  . ds  0  us bar  wr'ich  rapl  les  a  tendency  ‘  ■.'warxia 

<  *;ri vec f  1  on  »r  ’he  masne‘ic  Me  d  with  the  p.aama.* 
Such  t -jfivmct  ton  ‘emts  to  produce  string 

•  oncant ra‘ 1  >ns  of  current  at  the  entrance  and 
exit,  due  to  the  presence  of  a  low  bac* 


•l*ctro«otiv«  fore*  '.emf'.  From  Ob*1*  law 


J  *  JxjJ  -  o-»  *  <Jx8' 

tha  plaaaa  »u*t  taka  tba  nntlrn  alectric  ft* Id  in 
t.b«  Ab**nc«  of  back  emf,  thu*  ‘'orcin<  strong 
current  concent rmt lona  to  occur.  8eoau»e  of  t*u«, 
local  neatinj  of  the  ineulator  and  elnc'rod'** 
occurs,  accounting  for  their  aaaociated  ablation. 

"he  aia  of  true  paper  la  to  present  In'  a  •’-om 
an  experiment  Saaed  on  a  theory  'Pat  clause  the 
ability  to  ,-ontro.  '  he  du*riPut  '.in  of  'he 
magnetic  '’’.eld,  -,ence  urrent  denai'y,  m  an  MV" 
ar.-.et  channel  throuKh  careful  ax  .a.  varia'ion  jf 
tne  mtereleetroce  separation. 

The  Experimental  Apparatus 

A  joint  experiment  between  *  IT  and  95A 
Waehmgton  Research  Japoratorv  in  Alexandria, 
Virginia  wa*  conducted  at  'he  RDa  faci’.i'.es.  A 
•ICO  xJ  capacitor  lurk  wae  modified  into  -m  eigh' 
stage  20  ladder  ne'worX  to  provide  a  r_at 
'juaa  i -s  teady  current  pu.se  'o  tne  v<on  arcjet.  "he 
bans  hne  a  total  capacitance  of  J . mF  wnere  »hci 
section  is  separated  by  a  1'J  ad  solenoid  inductor. 
The  pulse  forming  network  °r\  19  designed  with  a 
high  voi'age  lgnitron  switch  at  each  stage  so  *bat 
the  output  waveform  can  be  modified  depending  on 
the  sequencing  of  lgrutrons  (refer  to  '’igijre  2-. 
If  a  single  lgmtron  ,  lgnitron  *8 !  is  fired,  the 
resulting  waveform  ha*  a  duration  of  about  0.8  ma . 
This  system  require*  a  matched  .oad  of  0...  0 
providing  a  maximum  safe  current  pul*e  of  J2  xa  to 
the  fFO  arcjet.  If  lgnitron*  *1  end  ex  are  fired 
a  i mu  1 1 aneou*  1  v ,  the  impedance  of  t-t  bank  drop*  m 
hal4'  and  a  pulse  of  about  J.4S  m*  is  achieved, 
•t-'b  this  conf  igurat  ion ,  a  maximum  safe  cu"tnt 
level  of  65  kA  can  be  supplied  to  'he  arcjet. 

Tbe  i*rV  i*  connected  *0  the  arcjet,  wmch  1* 
situated  ,n  a  cylindrical  stainless  steel  vacuum 
tanx  6  ■  long  and  0.6  m  m  diameter.  "he  vacuum 
i*  maintained  at  a  static  level  of  <  10  '  torr  bv 
two  six  inch  diffusion  puma .  '"he  'ank  and 
associated  pump*  are  electrically  '’’.oatiig  >0 
prevent  spurious  alternate  current  path*  'or  tb« 
-'fV-srcjet  circuit.  The  kvj  arcjet  aa*ejab'.v 
.‘rawing  in  figure  d  shows  the  relation  >■'  its 
various  component*.  "he  exterior  of  the  arcjet  -,s 
constructed  as  an  anode  sleeve  which  is  isoiate.; 
from  the  vacuum  tank  bv  a  t’lej<igla*  f'.ange.  "he 
sleeve  and  "ne  outer  portion  of  tie  anode  are 
hardcoat  anodized  wit.n  a  7  oil  layer  to  inhibit 
current  attachment  bevond  the  anode  lip.  "he 
cathode  and  an**  injection  assemblies  both  slide 
into  the  anode  sleeve.  A  boron  nitride  insulator, 
seated  in  a  “lexigla*  assembly,  isolates  'he  anode 
‘'r*  the  cathode  while  a  Wy'.nr  ''1.0  and  'ape 
jacket  surround  'he  cathode  *leev»  comp. "ting  tie 
anode/ cathode  insulation.  by  constructing  tie 
• eed*  to  the  arcjet  in  thia  coaxial  manner,  ‘he 
parasitic  inductance  of  the  device  is  xep  .  to  a 
minimum. 

“a**  is  provided  to  the  channel  'rom  a  large 
reservei'  placed  close  to  the  HL,n  er'-jet  hanne.. 
iix  valve*  feed  a  proexpans ion  ■  hamper  direct  .v 
up* tree*  of  48  cnoxed  orifices.  'he  gas  ' .owa 
* rough  the  choxed  orifice*  and  expands  through  a 
boron  ni'ride  insulator,  which  a. so  baa  Id  larger 
diameter  hole*,  into  the  channel. 

"he  anod*  la  mad*  of  aluminum,  primarily  'or 
it*  mech  inab  1 1 1  ’  y ,  low  weigh',  and  -nt  .  o*t  .  "he 
cathode*  ar«  mmde  'row  copper  pfiiari.v  •  >r  ’s 
e.ectrical  propertie*,  and  it  .*  both  lea* 


expen*ive  and  euier  to  machine  into  the  required 
geomet  riea  than  uther  *or<9  dasirabia  Mttria.*. 

*  ’’ektrornx  da'  a  acquisition  syataai  vs 
utilized  'or  la'.a  col  .action,  Mr  1  pul  at  ton ,  and 
'.i,  r»  )«  kept  inside  s  ergs  farsday  csj  e 

*  or  shielding  against  stray  e  L  rc  1 1  omegne t  i c  10  :■«. 
A  b  lock  d  l  ag  r  a*  of  r  h#  MHO  thruater  i y ■  t  as  vs  i9f n 
.n  'ig\jre  4 

'ruf  i*t  Ion  of  th*  discharge  occurs  after  *  he 
bank  is  «  her  god  fo  '  he  appropr  iat«  -ly^r  sting 
vh.jc.  At  '  tAOin*;,  the  valves  are  .jp^n^U  wi'm 
a  )  kV  pula*  an*1  are  ha.>l  open  on  t  l  .  tne  required 
state  anas  * '  1  j«  rate  .  s  r cached.  After 
’Ti*  dr.  av  iD^roxiMt«*lv’  jJ  sis  .  the  requ  i  r**d 
ire  *  several  <iiovol's  appear 

s*  ri)ii  ’hr  e**ct'0<ie  gap.  h  reaxd<.>wr>  ensues,  ino, 
*,rMr  a  abort  transient  period  f vpica.lv  .  hO  *s  , 
-UH.sv-stRat.lv  oorratv on  is  reached. 

Tnr  .later  vo.*a><R  aeasurrd  by  a  1000-  \ 

Trxtronix  voltage  probe  opt '  •  a . 1 y  isolated  *  rom 
* ne  ta^a  ac^u is v;on  svs’ea.  K.oarin<  po^rntia. 

^  Ti»*‘csufRd  rr.ative  *j  r  be  anoor  * :  os  s  si  sp.e 
'..ini^niu  i  r  probe  tbro-u^n  a  w1.1*'  1  vo  *t  iit(e  U  v  ;  :rr 
”ne  ♦".oat’nn  potential  probe  vs  constructed 
pr  vaanw  from  a  snort  piece  of  fhn  'un^sten 
The  total  current  is  arasured  by  a  .ir^e 
^o^owski  loop  surrounding  the  cathode  s.eevr  tt 
the  base  of  the  power  connection.  **b '  s  probe 
ueHsurm  the  time  rate  of  change  of  the  f*i'c:oi»*d 
magnetic  flux,  whicn  vs  integrated  to  giv*  t ne 
magnetic  field,  and,  from  Aapere  s  '.aw,  'he  total 
enclosed  current.  An  active  integrator  with  a 
J'7  us  integration  time  constant  and  a  70  mm 
associated  'roop  tise  IS  use<l  wrh  this  WogowSkl 
’■’*>p.  Knciosso  current,  at  points  inside  the 
is  <wmu>  *«1  h-s  a  see  1  f  0  turn  ^oiqmiu 
Ooo  enclose**  ,  >«  a  !ong  K’yrev  tube.  ^ajn,  sn 

a<-'  iv*  m*  sirs*  or  with  a  J .  .9  a*  integration  tis»e 

*  •instant  an*:  ^  ‘0  as  aeaociatsd  roop  ine  .s  usf  l 
to  m*rgra*r  *  he  1 1 M  derivative  of  the  magnet 
induction  f  e.d. 

■'or  »hia  experiment ,  'he  bank  waa  fired  with 
#!  and  ^8  vgnttrons  so  tnat  an  operating  point  of 
>0  kA  at  4 .  )  it/ 9  of’  Argon  was  obtained. 

"wo  wi*n  arcjet  «'nanne.a  were  constructed  and 
diagnosed.  A  constant  area  :hannri  cons  is  ts  of  a 
cv  .  ndf  leal  an*  Mia  and  concent  r  ;C  ca  *  h<K‘.e  with  a 
Cons*  ant  «i t •»  e  L  r-  '  r<  -  de  srpar a  '  i  5n  o'  ^  at .  s 

se.-..oc.  chauna  .  :.ne  eod :  ♦  ;  ed  channei,  was 

h/taed  on  a  .alcu.at.on  '*oa  ‘hr  *  heorv  v  h 
-"•'erencr  \  for  mfoni  current  drnaitv 
d. is' r ibuf  lona  along  the  r. ret:  odes.  To  ainimi/.e 
cv  indncal  »,7'c  s,  ’hr  mnular  channrl  geometry 
wiki  oasrd  on  a  .  arge  cathode  diiuaetrr  compared  to 

*  hr  r  tr  re  Le*:t  »*o«tr  gap.  *or  a  mass  flow  rate  of 
4.  J  g/ *  of  Argon  "i'  an  operating  point  of  about 
^0  ’xa,  s  chan  nr’,  wi*h  the  arrs  variation  shown  m 
^UTirr  ^  is  prr**  icted  bv  the  *  heorv.  In  order  to 
Rxaa’ne  'hr  hehavv  or  0*'  '  hr  *wo  channels,  we  chose 
to  operate  • n  *  he  vicinity  of  onset,  where  onset 
.s  predicted  by  'he  'heorv  as  the  solution  where 
zero  drnsL'v  vs  enroun'er**d  scasewhere  along  'he 
anoda.  "he  'Vsign  'heorv,  'or  which  this 
calfulatlon  vs  'ased,  neglect  a  'he  ini'va. 
■■’*ii?  a'  ion  nro«*«sses  'ha'  occur  s  reg  ■  on  of  the 
discharge  dmnna'rd  bv  "osic  heating  where  arge 
-  h'vngra  in  plasma  velocity  and  'ensvtv  occur.* 
'-'or  his  reason,  s  low  speed  near  the  entrance  *or 
a  given  mass  riuw  rate  earls  to  large  var  iations 

n  ‘he  l  n r  ere  te<  *  rode  separa'  ’.on .  This  vs 

vtBMd  lately  seen  in  rvgure  5  at  the  channel 
■  »■ -an*  a  ipe*>-a»  ^*osi  the  point  of  avniius  area. 
'(<»wm*vRr,  f  -r  "m  ex  peM  merit  ,  in  ordsr  to  ansur  • 
14m*  v  >n  at  the  inlet,  only  the  portion 
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downstream  fro*  the  point  of  minimum  area  wa» 
flared.  "''he  up*trea*  portion  remain*  'ixed  at 
. . H  cm.  Hence  thv*  1ev  .ce  i*  called  the  modified 
flarec  channel.  In  both  geometries,  the  anode 
remain*  the  sium  and  all  charnel  variation  .» 
taxen  by  the  cathode.  Thi*  l*  done  for  simplicity 
in  the  machining  of  the  electrode*. 

Results  and  Queue* ion 

Roth  channel*  were  probed  to  determine 
currant  and  ♦loatmg  potential  distributions. 
Typical  probe  ‘"eaponae*  are  ieen  in  figure*  t>  arul 
** .  "he  enclosed  current  and  floating  potential 
contours  tor  *ne  two  channels  are  founa  n  'igures 
R  and  9  respectively.  Vet  emmet  ion  of  the  plasma 
potent’,  a  1  f*rom  floating  potential  measurements 
necessitates  .Knowledge  or’  the  electron  temperature 
and  and  ionic  species.'  Other  work  has  shown  'hat 
'his  temperature  variation  is  saa..  so  that  'he 
differences  m  '1 oat  mg  potent is  sufficient 
for  the  ca.culation  of  e.ectric  ^i-lda. 

In  both  channels,  la^ge  ax:«.  urren's  are 
found  in  the  anode  region.  "his  up.  .*»a  large 
.ocal  values  of  the  Hall  parameter.  from  *.  v 
floating  potential  contours,  the--  a. so  appears  to 
be  a  *arge  positive  potential  jroo  at  both  anodes. 
\  .a rge  Hall  parame'.er  implies  a  a  deficiency  >f 
charge  carriers  which  a)rr»«.at*<i  w  th  ’he  ar«*» 
positive  anode  fall  required  to  pass  a.,  of  * \e 
current.  Based  on  reference  it  s  apparent 

'hat  both  channels  are  operating  beyond  the  onset 
lev-1. 

Typical  terminal  voltage  trace*  do  not 
display  ihe  high  ''requenev  oscillations  that 
•ypica*ly  characterize  onset.  Megahertz  voltage 
oscillation*  are  seen,  however,  m  the  floating 
potential  characteristic*.  "ermine  1  vo.tage 

scaling  shown  m  f,.gure  10  does  not  display  t.ie 
cubic  dependence  of  vo*tage  that  one  would  ••kpect 
to  precede  onset.1  •• 

"urchi  suggests  ♦'hat  an  increase  m  ablation 
as  the  exhaust  speed  exceeds  the  A.'ven  c-it  icsi 
speed.  baaed  >n  ablation  energy,  n«v  preclude  a 
cubic  dependence  of  volt eg*  o-n  cur  rent  in 
arejet  operation,  even  /  ne  dominant 
acceleration  mechanism  is  -  lect  romag^e  t  ic.  * 

Examination  of  both  thrust-’*  snow*  strong 
evidence  of  ablation.  Both  anodes  show  erosion 
primarily  at  their  exit  lip.  Both  cat  node*  .li.so 
show  erosion f  "spec tally  \n  the  entrance  region. 

n  ♦‘act,  after  lea*  than  100  discharges,  Copper 
f rom  the  cathode  is  plated  on  'he  boron  nitride 
nsulatmg  mjec'or  plate.  "he  anode  wear  i.a 
similar  .n  both  h»*'r,.est  however,  the  cathode 
wear  is  marked. v  d*. ff-rent.  for  the  constant  ar-a 
channel,  bands  about  a  ha.f  a  centimeter  wide  are 
seen  at  both  the  entrance  and  t  ie  -xt*.  "he  banc 
at  the  entrance  looks  like  *• -sb.lv  machined 
•■upper,  indicative  of  strong  and  even  -roaion. 

"he  band  at  the  exit  is  well  defined  but  is  not 
quit*  a*  clean  a*  the  one  a'  'he  entrance.  The 
''and  at  the  entrance  ar-  .srg-.v  composed  o*‘ 
semicircular  clean  regions  adjacent  to  the  too** 
. njection  site*.  *ore  c~etu j  sr'ots  are  noted  about 
a  half  to  one  centimeter  ' ur ' hr  r  downstream 
directly  m  line  with  the  mass  inj-ction  hoi-s  in 
'he  boron  nitride  insulator.  "he  remaining 
cathode  show*  strong  evidence  of  '*ine  1  me  nr*~ 
structure*  extending  'rom  the  cathode  downstream 
-nd  almost  to  the  cathode  root.  "h— se  fine 
structure*  are  branch- like  and  tend  to  branch  m 
an  up*tream  direction. 

In  the  modified  flared  channel,  the  clean 


band  at  the  exit  r*  absent,  but  there  i*  a  wid* 
clean  band  the  entrance  between  0.7B  cm  'o 

1 .  .1  cm  m  •xtent,  ^nlixe  the  or  her  bands,  'hi* 
<>ne  has  a  very  poor  ’  y  defined  downstream  wig*, 
"he  fine  l  t  rm  •truc'ure*  are  completely  absent  un 
' h : •  cathode.  Drawing*  of  both  cathode*  are  found 
in  figure*  Li  and  111 . 

>  igure  13  show*  a  graph  of  the  current 
density  along  both  cathode  surface*.  In  bo'h 
cases,  strong  concentrations  are  found  in  the 
initial  portion  of  the  channel.  In  'he  con*tan* 
area  channel,  the  current  density  drop*  'o  «  -  a* 
value  at  sid-channel  and  then  m<  rra*#»* 
dramatically  near  the  exit.  The  current  lens i ’ v 
m  r  he  modified  flared  channel  remains  high  x,,r  : 
j ua t  before  the  beginning  of  the  flar— 1  region 
where  it  subsequently  drop*  » ign i f leant . v  be.ow 
'he  level  found  in  the  constant  ires  channel. 

"he  theory  tenets  to  support  'he  trend*  seen 
n  figure  13.  At  the  channe.  erf  ranee,  the  piasba 
.3  do v  i n g  s  lowly  so  'he  back  em r  is  low.  . n  t •  ■  -  b 

Charnel.*,  the  ’.  n  fere  lect  rode  distance  at  *  be 
►"'trance  is  tie  same,  and  both  plasmas  must  t  i*e 
'  he  *otal  e.ectric  field  in  'he  .to sene e  <>♦*  my 
si  gn  i '“leant  back  emf.  At  the  exits  'he  ^agne*  . 

*  leid  ir  opa  toward*  zero  and  jnce  again  ’be  *’m<  « 

*mm?  i*  *mai..  In  this  region,  'he  p.asma  jnc- 
again  aua  ’  '  ake  'he  'utai  electric  *'ie.d.  h«»wever, 

*  »  •  •«#  i  f  ■  -■!  f '  s r*«4  -  hann*  i ,  t  h* 

'  w  tatfii  du*  to  the  im  rinsed 

i  n  t  er*  Lecf  r  o«ls  separation,  *o  *  a' rung  cur? -at 
c jncent r*t ion  i*  mitigated.  Since  the  current 
density  i*  luw*r  in  the  fl*r*d  region,  and  the 
*ajwe  currant  i*  pa**ed  a*  in  the  constant  area 
•bonne l,  'be  current  density  in  the  upstream 
porf  lun  of  the  modified  flared  channel  must, 
increase.  This  may  account  for  higher  ohsic 
heating  which  m*v  lead  'o  an  elevated  please 

*  -taper  ntur  * ,  and  'berefore  higher  iomza'ion 

act  ion  i/id  electrical  conductivity,  m  'he 
modified  f  «r  ad  channel. 

rhe  ♦ heory  also  indicate*  that  high  magnetic 
Reynolds  nuxOcr  flow*  in  constant  area  ob.inne.a 
♦end  'o  have  strong  current  concentration*  at  'he 
- v i '  as  the  magnetic  field  goes  to  zero.  Also  it 
snow*  that  'he  magnetic  field  in  a  r ared  channel 
s  sure  even  v  distributed  due  to  the  controlled 
p  asma  expansion  and  the  convection  o e  the 

m/vgn-' ic  f'.e.  d  with  the  fluid  at  h.g.her  magnet  ic 
'■rno  <t»  numbers.  In  this  opera'  mg  •“egime,  *  he 
magnetized  plasma  behaves  ona.ugou.3lv  ' o  an 
O'cinarv  coorpres*  ib  le  gas  except  that  'he  sp^ed  of 
iivjvi  is  replaced  by  the  magnet  oacoua  t  ic  speed, 
defined  as 


which  is  a  r<x*h  t n*f  ion  of  the  speed  of  sound  and 
t  ie  A 1  *"ven  m  . 

Data  'o'  'he  constant  area  channel  doe*  not 
seem  to  exhibit  this  behavior.  rh*  'neorv 
neg  ecf*  the  e!'*ect  of  viscosity  which  say  ac  «".uit 
for  some  of  the  difference*.  In  '.his  channel,  the 
n'-esence  o*’  a  large  axial  component  >f  the  current 
act*  'o  dev’l-ts  the  anode  and  presaurlres  'he 
-  at  hod*.  "h*  movmg  piasms,  blocked  bv  'he 
cat  nod*  *us^  *ove  along  it,  which  may  b* 

sj*'ici«nt  fur  viscous  dissipation  to  a  1 uw  the 
P  asm*.  f  the  p.asma  a peed  is  sufficiently 

low- r*d,  the  correspond  mg  magnetic  Heyno  vis 
nuimiar  for  t. h*  arejet  may  be  low  enough  for  a 
un^orm  axial  decrease  of  magnetic  f  .eld.  If  *he 
•negnetic  fi*’d  is  sufficiently  low  at  the  exit, 
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ordinary  gas  dynamic  choking  can  occur  a*  in  the 
cue  of  ordinary  ga m  dynanci  with  fanno  f-ow  for 
flowm  with  friction  and  Rayleigh  flow  for  ^low* 
with  Heatinf).  fhua  limiting  the  effective  exhaust 
veiocity  to  that  cor re* ponding  to  Mach  1. 

However,  the  modified  flare-:  channel  does 


exhibit 


of  the  trends  anticipated  by  the 


theory.  Although  pressur izat ion  of  the  cathode  is 
evidenced  by  the  current  contour*,  the  Glaring  of 
the  cathode  tenda  to  reduce  the  obstruction  of  the 
plasma  motion.  Also,  an  increased  pia aaa 
temperature  can  lead  to  ei  strong  drop  in  plasma 
viscosity  as  ihown  in  *fijfur«s  14.10  This  could  act 
to  sufficiently  lower  t*ie  viscous  losses  -x  pec  ted 
in  this  region.  There  for* ,  the  magnet,  ic  Reynolds 
number  aav  not  be  reduced  significantly;  m  rac', 


the  therroodvnanic  expansion  encountered 


in 


t  he 


nozz.e  ’-egioo  would  ai.ow  for  an  increase  m  the 
plasma  velocity  beyond  the  appropriate  sonic  or 
sagnetoacouat ic  speed,  which  would  translate  to  an 
increase  in  the  magnetic  Reynolds  nua^^r 

rhis  is  farther  evidenced  by  'he  difference 
m  the  terminal  voltages  measured  .n  '  ne  ar-.  _■*»». 
*at  *  mez-Sanchez  snows  that  the  i-'a;  v-x.'sge 
decrease*  with  increasing  magnetic  Hevno  .  '• 
n^jMOer.*  In  this  experiment,  ‘he  constant  ir-i 
channel  is  seen  to  operate  at  a  mean  vol’age  ot 
‘*4.1  t  1.6  vo*.ts,  while  the  modified  .a  red 
channel  operates  at  70. 9  -  l.S  volts. 

The  plas*A  discharge  near  'he  anode  appears 
quite  different  ^rom  the  cat  node.  rhe  anode 
current  densities  and  voltage  drops  are  plotted  m 
figure  15.  **oth  geometries  exhibit  .arge 
potential  drops  indicative  of  starved  anode*.  for 
the  most  part,  the  magnitude  of  the  voltage  drop 
follows  the  magnitude  of  the  current  density.  The 
only  piece  where  there  is  a  deviation  this 
-'ehavior  is  at  the  entrance  to  the  modified  flared 


i  nanne 


Here,  high  current  density  is  associated 


With 

the 


low  voLtage  drop.  This  tends  to  suggest 


deve lopment 


concentrated 


suostant  vated  through  examination  the  anode. 

In  'his  initial  region,  there  is  a  high 
•iesand  to  pass  current  due  to  the  .ow  back  -«*  , 
but  the  gas  has  a  low  conductivity  since  it  h/ia 
not  yet  become  fully  ionized.  rhe  .oca  l  -lee trie 
'’.elds  may  not  yet  be  sufficient  to  support  a 

diffuse  discharge  so  the  plasma  mav  have  no  other 
choice  but  to  arc  strongly  and  ocn, ,v  heat  t.ne 
p-.  asms.  as  the  p<asma  is  converted  -towns t  ream, 
its  conductivity  increases,  and  tne  arc  anv  se-owe 
more  diffuse.  Eventually,  no  severe  arcing  *.* 
necessary,  as  the  conductivity  and  "he  'reck  emf 
have  risen  to  a  level  so  that  a  more  diffuse 

discharge  may  be  supported.  This  is  evidenced  dv 
'he  wear  pattern  on  the  ms.de  of'  the  anode. 
Arcing,  adjacent  to  the  injector*  at  '  he  anode,  .  s 
seen  with  bracks  that  appear  'o  v>t  inu* 

downstream.  As  the  arcs  continue  downstream,  tw 
appear  to  broaden  and  weaxen  urfil  ' he  wwar  on  '  y» 
anode  surface  oecomes  more  azitrnt  ha  llv  'irdon. 
'his  occurs  about  a  third  or  the  wsv  down  the 
anode  from  tne  c nanne1  entrance. 

Significant  erosion  can  occur  m  other  ar*a* 
of  the  MP9  arejet  'or  reasons  "hat  are  not 
necessarily  associated  wi'h  an'.''".-  current 

conduction  onset.  In  this  experiment  erosion  o' 
'he  cathode  has  been  local  Iv  -.ler^-aaed  bv  ''Taring 
'he  cathode.  rheorv  suggests  that  ‘his  :*  also 
valid  for  'he  entrance  region.*  *v  d-c* easing  'he 
local  electric  field,  'he  curr^*  *  *P r-*^  ou' . 
thereby  decreasing  its  density  and  'V  associated 
oral  onmic  heating.  ’’hi.s  ^o*  on.v  pr  ••—.  '*  'me 
cathode,  but  the  insulator  region  as  w*  _  . 
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Preliminary  examination  of  data  'rom  an  MP 0 
arejet  experiment  «,nown  that  a^ejat  operation 
an  be  modified  through  a  variat.on  of  'he 
i  r  t  ere  ect  rode  separation.  '■*urf  h*r 

experimentation  should  sa  done  to  study  the  effect 
of  increasing  the  m t ere i ect rode  distance  af  'he 
entrance  m  a  similar  Banner  as  was  done  at  the 
exit  o  alleviate  the  high  one*” t rat  ton  or* 

current  a.ong  the  cathode.  rf  ’.he  expected 

decrease  m  current  density  is  -ncoun  t  ere-1 , 

additional  *xpsr imen t s t i on  using  be' ter  cat  node 
materials  with  superior  ablation  proper' \ea  such 
c*  '  hor  iat  e«i  tungs'en  or  barium  ox  i  he  .aprcgna'eu 
'ungs'en  s.hou.J  >  lone,  as  .  .  'is  :on*  .  -'uod 

a'ucv  o *'  the  'rfV9  -trejet  plasma  p--opwr''. es  so  .as  to 
md-'-a * iwicl  oss  merhanistns  such  as  'rict.on  in  'he 
.ire  jet.  And  ina.lv,  .uMltional 

-  x  per  i  men  t  ».t»  i  on  should  oe  conducted  to  s'udv  '  ‘>e 
*  -  r  i  n  oe  4  ow  onset  . 
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